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ABSTRACT 


A new approach is given by the author in this thesis to solve 
the problem of time delay and time of arrival in quantum mechanics. It 
is shown that time of arrival is an improper and not a proper quantum 
mechanical observable, since the associated operator is not self-adjoint 
but only Hermitian. This means that the fundamental property of self- 
adjoint operators, which supplies a complete orthonormal set of eigenfunc- 
tions with real eigenvalues does not apply. All the familiar results of 
quantum mechanics depend on this. Some authors have attempted to quantize 
certain functions of time of arrival in the hope of achieving self-adjoint- 
ness. A different method has been selected here. Instead of attempting 
to construct a self-adjoint operator, the author outlined a detailed proof 
of all orthogonality, linearity and transformation relations directly in 
the case of the Hermitian time of arrival operator without reference to 
the general results for self-adjoint operators. In this case the time of 
arrival operator is linear on subspaces of the configuration Hilbert space 


known as time domains. 


Also originated by the author, is the unique method for measuring 
time of arrival and hence time delay at a specified point in space, instead 
of time delay associated with some angular momentum quantum number 2& . 

The concept of eigenfunction point phasing is introduced and is justified 
by examples throughout the chapters of the thesis. Time delay is discussed 
for one dimensional transition and simple three dimensional radially symmet- 


ric elastic scattering. As the scattering differential cross section 
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determines the absolute value of the scattering amplitude, time delay is 
shown to determine its phase. Finally, some approximation methods are 


discussed. 


The introductory chapters prepare the way for the main theory 
to follow. Chapter I reviews wave mechanics and proves the uncertainty 
relation, and uses this result to obtain the minimum uncertainty wave 
packet. Chapters II and III discuss classical time delay in one dimension 
and three dimensions respectively. Chapter IV briefly examines the classi- 
cal conjugate Hamiltonian formulation of the time function. Chapter V 


begins the quantum theory of time measurement. 
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CHAPTER I 


NON-STATIONARY STATES OF A FREE PARTICLE 


A stationary state for a particle in quantum mechanics is 
an energy eigenstate with time dependence of the form eee jo teil rel 


Hae ee Deen Raat. : 
state satisfies the condition that [| is independent of time. For the 


one dimensional free particle equation 


2h ee 


Siete Set 
TC aba Dae? oe 
ox 
the stationary state solutions are Bees) where the energy eigenvalue 
nk? 
is E = hw = ayer The general solution to (1) which is nonstationary, can 


be evaluated as a linear combination of stationary states as 


Wee = + | A Gee dice: (2) 
V2 * 00 
A(k) determines the distribution of momentum components in the wavefunction 


w . The factor a is included for Fourier transform convenience. If we 
V 27 


define g(k,t) by 


5 ink? 
Ce awaChy ce oon (3) 


then > and wy are Fourier transforms of one another, that is 


W(x,t) = | h(E Wee vied wee 
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yw and are respectively known as the position and momentum 
wavefunctions. They are always related by the Fourier transform equations 
above, even in the case of a non-zero potential. However the time depen- 
dence of the momentum wavefunction does not have the simple form of (3) in 


this case. For a free particle we can write (3) as 


cae 
-ihk (t-t 9) 
2m 


$(k,t) = $(k,tp) ex 


which gives the time dependence of © at any future or past time in terms 
of the value of > at to - Since the Fourier transform is unitary, we 
can express wW(x,t) in terms of the x dependence of Ww at tg - Using 
(3) we can express A(k) in terms of o(k,ty) and hence in terms of 
W(x,t,) using the inverse Fourier transform. Substituting in (2) and 


integrating with respect to k we obtain — 


=, 2 
Pe adhat ™m = im(x-x) = 
WGet) = 90 /whCe=tp) [ hie DlaterMaace saat Cae Me se 


These coefficients come from the internal integration in the expression 
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ik(x-x) - (t-t,) : 
w(x,t) = 7 fs v(x, t,) if e ae ° ona B 


Observe that if A(k) and hence g(k) is square integrable, 


then so is wp(x,t) and moreover 
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[ |v(x,t) |“dx = | V (x,t)W(x,t)dx = | |o(k,t) [Zak = | Jack) |? ak . 
—0O —o —OO —0O 
If all these integrals equal 1 , the wave packet representing the parti- 
cle is said to be normalized. To show the norm preserving property of the 


Fourier transform which establishes the above result, we substitute the 
. *° * e * * 
Fourier integral for wW in terms of > in the expression ww dx 
‘ —0O 


using two different integration variables k and b inthe 6 function 


integration. We can then separate out the delta function integral 
1 [ _ix(k-b 
6 (k-b) Yeu | aes =>) dx .yerine: Symbol 6 te not a function at all, 
—0O 


but 6(k-b) is the kernel of an operator (the identity operator) on the 
linear space of Fourier transformable functions. For this we see that the 
Fourier transform operator F preserves norm. Thus the normalization 
integral is independent of time for both position and momentum wavefunc- 
tions. 

kx 


Using (4) with a 0 if we substitute w(x,0) = Be and 


ikx - iwt ; : ‘ 4 
integrate we get wW(x,t) =e - which is the travelling wave which 


is not square integrable. On the other hand, if we take the wave packet 


localized at t = a 0 , namely wW(x,0) = 5(x) we find 
(x,t) = 4) ane exp CE) which is not normalized since 


[ 16 (x) |? dx = 6(0) which is not finite. We would have to take 
—0CO 
W(x, 0)= Vv S(x) to get a normalized square integrable function, but this 


would give w(x,t) = 0. To understand the reason for this, let us look 
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the integral in (4) and evaluating (details omitted) we get 
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From this equation (5) we can see the way a wavefunction Spreads as a 

Eunetion of timeydepending on itssinitial configuration.Meeit ont >> m 


we have 


1/4 (ihot)-1/2 im’) 


ve (x,t) o (= ) = exp (Fe 


Thus for large a, the initial configuration (which is normalized) more 


closely resembles Vop(* 9) = ¥5(x) and the wavefunction W(t) 


approaches the form exp (=== 


TT as Q-7> © with amplitude degenerating to 


zero. This explains the reason we found (x,t) to be zero previously. 


It is apparent that the wavefunction by Gt) spreads with time, 
and has minimum spread at time t = 0. For large a (initially narrow 
distribution) the wave packet spreads quickly, while for small positive a 
near zero, it spreads slowly in time. This will be shown quantitatively 


in the following discussion. 


The expectation value of a given physical quantity A ina state 
ere 
wy is represented by <A>, = | w (z,t) A pCx,t) dx, where ) is normal— 
—00 


ized. Similarly, with respect to the momentum wavefunction 
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Vv 
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i @ (k,t) Bd (k,t) dk . A may be any operator on functions of x, 


and B any operator on functions of k . For the one dimensional free 


particle, some of the more important operators are listed below. 
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TABLE I 


OPERATOR REPRESENTATIONS 


Physical 
Quantity 


Position x 


Momentum / fh =o 


Energy = E 


Any of these operator equivalence can be verified. 
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| W(x) A W(x) dx = 


Position Momentum 
Representation Representation 
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The condition 


to be equivalent is 


fe "(k) B O(k) dk 


a ik 
for allefumetions=*o(k) “such that°e(x)=-—— o(k) eo * dk or 

27 °-© 

=] 
y = Fd where F is the Fourier transform operator. This means B = F AF 
since F in unitary. 
Now return to the function y (x,t) . The spread of a wave- 
a 
packet is equal to Hea oes = gee = eee Bue <x> = 0 “for Vy and 
Dee ae, 


so (ee = <x > = 


2 Ql im se sam 


This explicitly shows that t = 0 
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the time of minimum spread, and that narrow distributions at t = 0 


(large Q) spread more rapidly. 


For a free particle in one dimension we can demonstrate the 
following results, directly from the time dependence in (3) or (4), and 


the operators in Table I. 
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<s ile = ie nKy a constant independent of t . Moreover 


CO 
k. = <k>, = | LACK) |? dk , and may be positive or negative. Define 
—O0O 


P 
v= au . Then, 
m 
<x>, = v(t-to) + <x? 5 
0 
2 Z ae 
h 0 fhok 
Ee oe ero a ae a a tae where E, ? 0 9 
ox 
2 
iy ee 
and E. Saree aee gah s 7e Oi, Degee ald. Ey are positive real constants inde- 


pendent of time. E. is energy of translation, and the subscript does not 


refer to time. Ey is energy difference. 


If we define (Ax)? = <(x-<x>,)°>, as the spread of the wave 
packet, then (Ax)? is a positive quadratic in it achieving a minimum at 


some time which is known as the time of minimum spread. We can show 
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(Ax), = ae + as E(t t,) 


by using the time dependence equation for the free particle derived earlier. 
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Uncertainty Relation and Minimizing Wave Packet 


The spread (AA)? of an operator A on a given wave packet is 


: 2 
defined by (AA)” = <(A-<A>)*> . From this we can state an important 


theorem about Hermitian operators. 


Theorem: (Uncertainty Relations). 


Suppose A and B are Hermitian operators and [A,B] = AB-BA = 


iC. Then AA AB >F ce leh 


Proors (1) We first show 5 es = Im <AB> . Since AB-BA = iC we 


have <AB> -— <BA> = i<xC> . Also <AB> 


(W,ABY) = (AV,BY) = (BAY,Y) = 


* * 
(W,BAY) = <BA> . Hence <AB> - <BA> = 2iIm <AB> = i<C> proving (1). 


C25?) et Go = Al <A>seiand (66 @=uB = <B>). Then: a ‘and 1/6. are 
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Hermitian and moreover (AA)~ = <a > = <A> — <A> and 


(AB)* = <@7> = <B*> = <B>”. Then 


(AA)? (AB)? = <a><B2> = (yah) (W,8°V) = Con, om) (BY, BY) 


Low] |* |]ev}|? > |cap,ew|* . 


Hence AA AB > |(o,B¥)| = |(@,a8~)| = |<aB>| . Then 


(A-<A>) (B-<B>) = AB - <A>B - <B>A + <A><B> , 


ap 


and so 


<aB> <AB> — <A><B> - <A><B> + <A><B> = <AB> -— <A><B> . 
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|Im<AB>| , 

since ~A;5 (Hermitian imply “<A> , <B> real, Since 
iL a 

| im<aB>| = | <c>|= 5 |<c>| 


we obtain AA AB > = |<c>| : 


We observe that the equality AA AB = = |<c>| holds if and only 
if the following two conditions hold: 
(1) ap and §fW are scalar multiples of one another, 
(2) Re <AB> = <A> <B>. 


We can apply the theorem to the product Ax Ap. Considering 


the momentum wavefunction ¢(k) at fixed time, we have that the operators 


Prana, ew, ote : 3 
Seat) as aeand = hk are Hermitian, and ac =. ihi—— - 
3p 5k P [x,P]¢ ap SPO? 
pih = (>) = iff , so we have the operator equivalence [x,p] = ih. By 
the theorem, Ax Ap >a. To minimize the uncertainty product we require 


that (s-<x>)@ ard §(p-<p>)¢ be scalar multiples, and Re<xp> = <x><p> . 
Hence to minimize Ax Ap , o(k) must satisfy i ae ok), = <x>o(k) = 
y( (hk - Po) o(k) ] , for undetermined complex number y . We ignore the t 
dependence of ¢ ; Ap Ax spreads out to become larger at all future and 
past time and we assume t = ty is the time of minimum Co - solving 
the differential equation gives 

| B 


-i([<x> - YPg]k +.yh _ 
o(k) = Ne 9 


” 


pycnugnerl 
wn : 

Conti = 

pe ae 

7 oe 

(aol fe that ahaa 


vino boe 22 sbLoul <om| 3 = GA AS vobtaups add Jad opened ag | , 


iia anutothabs ona aebrotbet' oF 


qoujons sho Xo wolgislam wfase w7e. 44 hoe ~ Diy! 
7 ag 


. > €0> « Casi 


gnizeblano) . qh xf seubosg ada 02 eros of? tiqqs ond a 
s7o927He0 Of) tad3 eved aw ,Smid bowss 76 Gly eolssavievew ev2ae 
— toa) - ft = o{q,x] tne ,nalttnzat ats a4 <q han 4 t ~5m 
yé . a = [Ges] soneleviope a9zeyego oda send ew of , “— wt 6, 
extuper ow Joubowq elabbdtoom ada gulginta cT Ge 4 gh 2b osoed3 ® 
. 2Goe%> * <aeret bane golekatue tafuoe od p(éqe-q) ban (conn). 

« Gaigees - OE get yhetron sem CDG « qb xh oxlalate 09 | i 

2 si osongi ow. ¥Y sade zolqnos baniertetalmy wo? , Let = ar 
bao ovute} Tin te veyxzel omoned of juo ebagvye x*A ch ¢ § 20 — 
qaivioe . S(xh) Cumbnka Xo omta ay Bk 3 + 3 mania wet tale 23 
eovlg aotseups | 


for a normalization constant N_ such that [ lock) |? dik = oh. “inis 


implies Im y< 0. In addition, © must satisfy a number of other 


conditions namely 
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[ O°(k) tk O(k)dk = pp 


If we write yY=a- if where 8 >0O , then 6 becomes 


+ <Ok - py ak] - B B(k-ky)? 
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where Po = hk) 
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Applying the normalization condition we obtain N = ees . 
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otk? Bh 2 
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We now consider the final condition 
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Dx ip. yoinces Ax Ap Bis minimized at time t= ti > (6) represents 
o(k,t,) Using (3) we obtain 
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Equation (6) gives the momentum wave function at time t 


namely 
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The position wave function is determined by Fourier Transforma- 


tion to be 


V(x,t)) = W(x) = +} o(k) eM ak oor 
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The constant §$ > 0Q is related to the stationary energy Ey ‘ 
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en. 3) 
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wavefunction Ak Ax = 7 » we have Ax = ree Hence the parameter 8 


determines the stationary energy and the spread of position and momentum. 
We observe that as time t differs from th » Ak will remain constant at 
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Z 
(Ax)e = AB + (t-t,)* for B>o . (9) 
Equation (9) shows that for large 8 , the position wavefunction wy has 
a large spread and its time dependence is slow. For small 8 the wave- 


function wW is very narrow at th » but spreads out rapidly. 


We saw that the momentum wavefunction for the minimum uncertainty 


product packet was given by 
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O(k, 6) = ES exp[-i<x>, k - EB (k-ck>)? - os (t-t,)] 
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Here ty is the time of minimum spread Ax , <k> is the expected value 

of momentum and > san is the expected value of position at time ty : 
1 

8 > 0 is a real parameter designed to specify the relative amounts of 


position and momentum uncertainty. 
The position wavefunctions wW(x,t) can be obtained in two ways: 
(1) Fourier transform the expression for $(k,t) above by taking 
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and use the general time dependence integral 


eB Y, fhe Pes eee 
W(x,t) = Mises) eS ThCe=t) | v(x,t,) exp Geese? dx 


In either case, the result turns out to be 
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If we set <x> 202 ot. = 0, <k>se20. 248 aes we can reduce this to 
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equation (5), namely 
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This gives the complete minimum uncertainty wavefunction in one 
dimension for the free particle in terms of 
(1) t, the time of minimum spread, 
(2) <k> the mean (average) or expected momentum (constant in time), 
(3) <X>, the expected position at time ti > 


(4) Ak the momentum spread (a constant in time). 
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The position spread is (at time t) 
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The results in this section can easily be generalized to three 


dimensional cartesian coordinates. One dimension is used for economy of 


notion and simplicity. 
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CHAPTER IL 


TIME DELAY THEORY IN CLASSICAL MECHANICS 


Consider a rolling frictionless ball in one dimension, with a 
height component y and a position component x ina gravitational field 


g . Suppose the ball has to pass over a hill y = h(x) between x = -a 


ands x = tay. 


Figure 1 
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Assume the ball is a point mass and has no rotational kinetic energy. 
Assume that before interacting with the potential V(x) = mgh(x) the ball 


is approaching the potential from the left with velocity Vo where 


> mvs > MEY where Yo is the maximum of h(x) for all x. Suppose 
that the ball has position Xo at time to where Xo <-a. 
The x and y components of the ball satisfy y = h(x) . Con- 
; ; ee 2 ee 1 dx dy, 2 
servation of energy gives 7 MV9 = 9 ml (7) ce pe leer Meyesy, aso 
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In particular, for h(x) = 0 we get the free particle equation 
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ea 2 Ege eae) for the one dimensional classical motion. The 
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cates this is undisturbed motion. The time delay At =t-t" at the 
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point x is given by 
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(At), = t-t' = | 


Now consider the particle in one dimension which is under the 
influence of a potential V(x) . Unlike the gravitational problem above, 


this is strictly one dimensional. 
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Assume that) V(x) sis non-zero only for---a < x < a s» Let 


Xo < -a and assume the particle is at Xg at time to and moving to the 


right. Conservation of energy gives us 


2 Bee be dx 
7 MV) = > MX avi See eSOntnatl, ot l= to + | 


The time delay is 


(A =a f ee (2) 


The Three Parameter Potential 


In order to determine what information time delay measurements 
give about a potential function we introduce a model potential V(x) known 
as the three parameter poential. Results of measuring At aid us in 


obtaining the three parameters and consequently the form of the potential. 


V(x) 
<- 2|w- 
3 a 
> xX 
<d- 0 = (al Se 
Figure 2 


Piecewise Smooth Three Parameter Potential 
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V(x) is assumed symmetric about x = 0. The potential acts between -a 
and a where a=wtd. d,h and w are the three parameters needed 


to specify V(x) . 


Let us calculate the time delay for the three parameter potential 


in the’case of the particles under gravity. |Substituting hx) = Mo 
mg 
in the expression for (At), we get 
L + [V'(x)]*/m’p? a0) 
(At) = fe 
x BS eae ) Vo 
Yo TVG 
Lee us take x, = —a = -w-d and "x = a= wid. The integration over x 


0) 


can be split into 3 intervals. The eventual result is 


At is the time delay after a complete interaction with the three parameter 
potential. If we hold d and w fixed and calculate the limit of At 
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The time delay for the free particle in one dimension which inter- 


acts with the three parameter potential is 


at = HH - S42 E - / vi - 2 | (4) 
2 2h 0 
Ves =a 
m 
This is the same result we obtain by letting @->® in (3). It corresponds 


to increasing gravitational strength g and decreasing height h_ so as to 


maintain the same function V(x) for the potential. 
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The Smooth Three Parameter Potential 


—*+qd>-> ——< 4g >— 


Figure 3 


Smooth Three Parameter Potential 


The smooth three parameter potential is a symmetric potential 
function constructed from parabolic arcs. On the interval [w,a] on the 
x axis, V(x) is the parabola passing through the points (a,0) and 
(w,h) in the x-y plane which is tangent to the x axis at x=a. 
Similarly V(x) is defined on [-a,-w] . On the interval (-w,w) V(x) 
is defined as a parabola that makes V(x) and its derivative continuous at 
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The explicit expression for V(x) on each interval is 


V(x) = (x +wet ane on [-w-d,-w] , 
VGoe= i+ m = x” on [-w,w] ; (5) 
V(x) = 2 (xw-a)” on [w,wtd ] : 
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Let us determine the time delay At for the one dimensional 


interaction with the smooth three parameter potential. It is given by 


+d 
ae [ dx _ 2¢wtd) 


-w-d f2 2 Yo 
Vans V(x) 
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We assume Vo is large enough so that the argument of the square root is 


positive. The integral on the range -w to w is symmetric about x = 0 


and so 
+d W 
be = 2 | Rieter IL ERE ee oe  - S. 
W hee aed i AID eae ae 0 oe ot 2 Wh, Seeliee es 0 
Yo a(x w-d) Vo neh + 7 q ] 
md 
2h 
a 2 (wtd 
Act = 2d peeaeee in 2 + 2 oe sinh : eee (6) 
“6 indy 2) aye, 0 
ZV i os 
One could curve fit equation (6) which describes At as a func- 
LOU (OLS OV for w > ee ab to experimentally measured values of At 
6) 0 m md 
in terms of v,. This would give a means for determining the parameters 


0 
Wows d. ania ae 


s2 s0o1 oxmuipe wis To vibes a atl pk onbAagr wR ge 


O « @ dood Sizsovmve ef ow 02 we Sggny sif7 no Tezgesnt sar 


Zi. 


High Energy Approximations 


: : 2 2h 2wh 
pe en 2 retake i 
For large Vo satisfying Vo ss + mq we can approximate 
in (6) above to get 
de ~ fat Sup +e (7) 
md Vo 3m Vo 


Thus for high energy particles At is proportional to = and the con- 
Vv 
0 
stant of proportionality, given in (7), is expressed in terms of the three 
parameters of the potential. This only gives us one equation for the three 
parameters, if we measure At for large Yo and determine the coefficient. 
To get all three parameters we need three equations, and consequently it 


is necessary to expand At to a higher order. As an example let us take 


the piecewise smooth 3-parameter potential whose time delay for interaction 


is given by equation (4). Expanding in powers of Vo we get 
Ss) 2 
ae w C2wtayn , Gwiayh? » OF FZ DRT 
= 3 2, y m2 vs 3 
mM Vo m VG 0 


By curve fitting this expression for large Vo we get three 


coefficients, which lets us solve for the three parameters w,d and h. 


In general, measurement of At as a function of Vo determines 
some of the properties of the potential V(x) , but not the complete func- 
tion V itself. The time delay is unchanged if the potential is displaced 


linearly along the x axis, or reflected through the y axis for instance. 
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Let us suppose that V(x) is non-zero only for -a<x<a. 


The time delay for interaction is 
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ar = | 25 ROC B 
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Vv 
[v' = = v(x) 0 
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We write the integrand as 


Ce = 1/2 ae 2V(x) 1/2 
(vo ay m V(x)) a Vv (1 aa D) ) 
0 mv 
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“0 n=1 Gm@iyint2- m ve 
1 ° (2n-1)! V(x) an 
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Let ae = rf [V(x) ]” dx . We get At expressed as 
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o0 (2n-1)! V 
At = ) ———__4_.__ .. 
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By measuring At as a function of Vo we can determine the coefficients 
V. which in turn give us information about V. The total information 
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which time delay measurement gives us about the potential is the set of 


all coefficients Whe 
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Hamiltonian Formulation 


Consider the point mass in one dimension approaching a potential 
barrier V(x) om <a ~ x * a, with sufficient energy to overcome it and 


pass through. The Hamiltonian function representing energy is 


2 
H = spre V(x) . The time function Ty (ts P>x) is T(t, P,x) = t+ T(p;x) 
OH O oH O 
where Cae = so =o =]. We take 
x Op pao 
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x 2 
ieee ane vest enGpy = -@1? | [B+ V(x) - Te eee 
x 
ab 
T)(t,P,x) =t+tT(p,x) . It is easy to see that 1 satisfies 
, Otero ir ; ‘ : ; ; ; 
VOC) co. ae 1 which is the conjugacy equation. (T is conjugate to 


H.) We measure time t as t= To - T(p,x) where T) is a constant. 


This gives us the time variation in terms of position x . The momentum 


p is always determined from x and is positive, so p = v 2m(H-V(x)) 
dope LUCE LOMUr OG ese FOr constant, Hi) Thus 6 = To ~- T(V2m(H-V(x)) , x) 
gives the time t as a function of position x in terms of constant 


parameters To) and Hi. Substituting explicitly for ~t(p,x) “we obtain 


Xx 
r= t+ pi? | [H= VEN] (* ae 
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In our original notation, Tn to A x = Xo and H = 7m ie - At time 


to the particle was at position Xo outside the potential range and had 


speed v approaching the potential. We had the equation 


0 


telonsj0q 6 anicionorqge sore 
bias sb scooveve 02 ygzanya 


” 


al ygton. grbaasaaud 


(eagdt + 9 = Cen 3),7 oat Gi ater wolsouwt eats sdT paw +s 


siod ov r~ Rs -a5 oem 


£ ‘ ; . 

AP. . : _ 

. 36 1G ev ‘ Ey a= ga. 
x 


YT 


buns 3B 1 = ta + “ij i SS) = = (aq)? Indd anne 
i" 


gatiaijsas +t Jeg coe OF YenS 2) dl, % (xe) + 2° 


76 
od sgeguimes 22 1) .mokynye yooyrteoo ody ei sotdw is a2 


(} 


idaayenes «at 6 asdtiw (Ge.q@)l =)? 1 a5 3 ae ann 
4 


i " 


= 8 


ri 


ot 


r| 


a 
or 


my oe) 

r 
{ , 
ny 


suinoewem edt . x noliladg 46 tasee et nolssivey onis si2 si eavig @ si) 


te , Uev-tpee i? — 1-48 2 aedl . Haneda 10t x Je coksoned & 


0 


qbdane: Yo 2teoi ak x molaieay lo solganvt os an 2 See md? % 


-o ‘ieido Ow oCeyu)y 182 “falallqeo guigevitiedse .. 0 bth a” 


» aw ‘| Ni tm 

‘ (> 
gat tk guage bat gt gk «igh * gt amet 
cua SgnGr. Intabyveg thy salad gt Aetstaoq. ae oat 
nok snups oti oid it Gintitotaq. ora gaia 


y YO ea 


(ie) 7-H) =q of ,svilleog al bon * wtal bontevtob ayewla at 


_ 


ater 


vo: 
Vi, : 


0 


- 
—— 


. 


d ‘ 
t = to + Ape a from which we obtained equation (2). 
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We summarize the three parameter potentials as follows: 


(a) Smooth three parameter potential. 


V(x) = 2 +wrt aye FOL = wad <x we, 
d 
wh | lee 4 
V@)3= h + Sry ae pes for ewes xe we, a = wid; 
h 2 
V(x) = az (x -w- d) tom 7 wos xe wid 


(b) Piecewise smooth three paramter potential. 


V(x) = 2 (x + d+ w) fox S-w-d) <x wee 
V(x) =h Ove ew . a = wtd 
V(x) = 2 (=x +d+w) for (wis x < wid 

(c) Piecewise constant three parameter potential. 
V(x) = b for —29 — x ~ ise ; 
V(x) =h for =S5 x 553 : aus 2 S5 Daun 
V(x) = b foc Sikes 


In each case (a), (b), (c) above, we can evaluate va Vy “ V5 
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Nes = | [V(x)]” dx to get the results 
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If we know At as a function of Vo for darge Yo Subriciently well to 


determine the three parameters Vy . Vy A V3 » then the form of the poten- 


tial can be determined according to (a), (b) or (c) above. 
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As an example, letsus “take case (b) and solve for w,h.,d 


in terms of vy ; Vv, ; V3 . We have 
V V V 
edt , S-fatw , 3-44 0% . 
h h 
on V5 a V5 
Eliminating d and w we get d = 2(— - —) = 3(—--—) so that 
h hn? h 2 
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vs h - 3V., ne- 2V., = 0 so we can solve this quadratic and get 


2 
BV any ON, gavin 


h = eer ae pe aes ea. Knowing h , w and d can be obtained as 
1 
V V V 
Bey oe aes ey ee 
d = Tees =e and w 5 (Gr dja. 


a 
Determining the Potential V(x) From the Integrals see = | [V(x) ] dx. 
-a 


We saw that time delay measurements in one dimension yielded the 
numbers ie as the total available information about the potential V(x) . 
We shall see that ae do not uniquely determine V(x) , but do uniquely 
determine another function G(h) from which we can get all possible V(x) 
with moments Ve . Conversely G(h) itself determines the a - Gh) 
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Figure 4 


Constructing G(h) 


Let dx = ldx, | oe | dx. | (both as positive 
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ve = | [Vix a> Sdsoe) a cannot be specified arbitrarily. For instance 
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Vow sole andelie Nixie 07, Vout] 2 0. In the equation 


> 0 and 


nh 


max 
V = | h” G(h) dh we observe that G(h) > 0 and so V5 
y = n 


min 
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Let us now suppose the potential V(x) is an even function of 


x (symmetry about x 0). Suppose also that the value of V(0) =a. 
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V(x) 
gives 2x = | GQyyudy. This holds af a = V and V(x) is mono- 
o max 


tone increasing on [-a,0] and monotone decreasing on [0,a] . This 
eivesmecueda terms (ota) Go) m= =rorm=xec™ [=a 0)... sinees Gi) 0) Vand 

V(x) 50°) By symmetry we obtain x in terms of VG) for xe [0,a} 
In this particular case we can construct the unique symmetric even func- 
tion, monotone increasing for negative x and decreasing for positive x 
in terms of the function G(y) . Observe that this is the form for all 
our three parameter potentials, and that it is uniquely determined by 


time delay measurements. 


Since V(x) is assumed to be of finite range V(x) = 0 for 
Sua 4), So a, awe see that. Vie) = 70) and Vv = V(O) for this special 
min max 
symmetric function. These are the consequences of assuming that V(x) 
has a maximum and is symmetric about it. The monotone assumptions on 
[-a,0] and [0,a] make a single integration of G(y) over y _ possible. 


This gives the basic reconstruction of V(x) for a positive potential. 
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Let us consider one very special example of a one dimensional 
potential for which the time delay can be used to determine some of its 


properties. 
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Figure 5 
Triangular Potential 
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three parameter potential whose time delay is given by (4). We set d=a, 


w= 0 and h=h to get 


dt = - 242 fy -Y ve - By . (8) 
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Now let us suppose we are given equation (8) describing At asa 


function of Yo? and we know nothing about the potential V(x) . Expanding 
At as an infinite series in Vy we get the following. 
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We compare this series expansion to the one for the general potential, 


namely 
oo (2n-1)! V 
Cte ) n-1 = OS z 
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and from this we conclude that V_ = V can be expressed as an 
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CHAPTER III 


SCATTERING FROM A CENTRAL POTENTIAL 


b = impact parameter 
ee = scattering angle 
6 = initial (and final) speed 


Figure 8 


Central Potential Scattering 


A three dimensional classical scattering problem with a potential 
V(r) depending only on r , the distance from the origin, can be reduced 
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The equations of motion are LS Cs = a and a3 es = = = 0. These 


include conservation of energy and angular momentum. The conjugate momenta 
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In computing time delay, we assume that hard sphere interaction 
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One way to solve this problem is to construct the most general 
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The conditions required for interaction with the potential is 


Pans <a Vins amples as a tecessary “condition “b= a . “Since the 
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For a positive potential we have a > Le ee >. by fitiucenerale eAt sthe limit 


the denominator of the last integrand becomes zero. The result is 
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i is less than Se and is the radius of closest approach 
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Figure 9 


Plane Integration 
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h(s) : 
The function -————— is to be integrated over the region R 
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shown on the previous page. One can cover this same region by performing 


the 8 integration first. We get 
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For u?A,a (u) = u and the time delay is zero (no intera- 


tion). This shows us that the constant of integration is zero and hence 
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means 0 (bs Vo) is negative in general and we have a time advance rather 
than a time delay. This is not surprizing since the free particle path 


is longer than the actual path. 
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The potential V(r) must be such that a(r) is monotone increasing in r , 
invertible, with differentiable inverse. If Vo and r are sufficiently 


large, this is true. 
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Impact Parameter - Scattering Intensities 


Suppose we have a uniform beam of particles incident upon a 
Scattering central potential. Suppose furthermore, the beam has intensity 
15 particles per unit area per unit time when directed on an area ortho- 
gonal to the beam. Suppose T.(@) particles per unit solid angle per 


unit time are detected at a scattering angle 9 from the initial direction. 


ag Ca) 
F : é ; 3 Ss 
Then the differential scattering cross section is -—~= where do 


dv Ip 


is a unit of area at a specified location in the incident beam, such that 


all particles passing through dO pass through the specified element of 


solid angle d located at the scattering angle 6. Because of radial 


oO 

dependence in the potential, we have no $ dependence in scattering. x 

as a function of 9® can be measured experimentally from 

T,(9) Alternatively, this can be used to relate the impact parameter 
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to the scattering angle, and gives us a means to get 0, (b) for. fixed Yo ° 
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than v. and b. However with this additional information, 9, can be 
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obtained as a function of b , and thus the results in (10)can be applied. 
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Total Reflection or Transmission 


Let us consider LOS, where b= 0 as a function of Vo - 


In this case we return to the problem in one dimension. Either the 
Parcicte is )rerlected back along its initial’ path 0, =" ,.Or it over— 
comes the potential and passes through 0. = 0. In either case 0 (05 v5) 
can be measured as a function of v, and this provides information about 


0 


the potential. Substituting b = 0 in an earlier expression for dt we 


get 
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dimensional problems. Let us assume reflection (r_,. 7 0) . Let us also 
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We change the variable of integration in this last integral. 
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For this to hold A must be a very special point uniquely determined from 


the potential. We require  V(r)s= 0" forg rc > A and V(r) > 0 £for 
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This solution is complete except for one difficulty. The range 
A of the potential must be known in order to complete the potential V(r). 
This same problem appeared in one dimensional classical time delay in 
scattering. The integral involving A can be separated out and directly 
integrated. However first we can integrate the entire equation over u_ to 


get 
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the constant of integration determined to be A from the condition u=0. 


Now integrating the coefficient of A we get cancellation of A and 
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Notice that A _ has cancelled out of the expression and so the difficulty 
of knowing the point at which the potential vanishes has disappeared. (11) 
can be treated as a general result for monotone decreasing V(r) that 
approach zero. Equations (10) and (11) also indicate the importance of the 


impact parameter time delay Os and the ease of inversion. 


Example of Time Delay and Potential 
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Thus we obtain the explicit time delay (which is negative, meaning a 


time advance) for this example as 


Note that v1 and c are parameters determined by the potential, and 


are part of the definition of V(r) . Now consider substituting this into 
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The Scattering Angle and Its Inversion 


The objective in this section is to determine the potential V(r) 
given the scattering angle 6, asa function of b , the impact parameter. 
As we have seen, this function can be obtained by experimental measurement 
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We assume that the particles all have the same mass in the beam being 


scattered, and all the same (vector) velocity with magnitude Vo ° Define 
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proves that the unique potential with the given scattering angle 0, (b) 
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CHAPTER IV 


TIME OF ARRIVAL IN CLASSICAL MECHANICS 


Consider the one dimensional Hamiltonian with the properties, 
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Hence the most general solution is 
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CHAPTER V 


SCATTERING THEORY AND THE TIME OPERATOR 


Consider the Schrodinger equation, 
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Normalization 
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Energy Degeneracy 


In solving the equation HY (x) = tu), (x) we may have a number 
(in some cases infinitely many) linearly independent solutions. For a 
given wavefunction (r,t) only one solution of this equation contributes 
to the superposition over w in W(r,t) . <A solution i) of 
HY = hwy’ which is orthogonal to De » the contributing solution, is also 
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(b) Energy and time amplitude functions. 
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Observe that the position or momentum wavefunction completely determines 

all knowledge of the particle, but this is by no means true for the energy 
(or time) amplitude functions. Specifying the position wavefunction deter- 
mines the energy amplitude function but not conversely. Also, observe that 
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the energy amplitude function may depend on To 
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Interpretation of the Time Amplitude Function 


A particle may be detected by placing a screen in a specified 
position. Consider, as an example, a free particle in one dimension with 


a 1ocaiized positive: momentum distribution (i.e. o(kjt) = 0 fer “kK < 0). 
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Figure 10 


Positive Momentum Disbtribution 


This particle is moving in the positive x direction. Suppose 
we place a screen at x = Xo > and we ask at what time does the particle 
hit the screen. We may make a series of measurements and find a distribu- 
tion of different times. The probability that the particle hits the screen 
at a time between t and ttdt is [rct) |? dt where T(t) is the time 


amplitude function derived from Xo phasing. 
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Let us write T(t Xp) for T(t) “with Xo phasing. We see 
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the free particle Hamiltonian. 
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Definition: A Time Domain D is a subspace of the Hilbert space spanned 
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Then we have 


-iw(t-t,) 
> if * — 
T(t,r 9) = f seh us * r)e 0 aw] Wr, ty) dr 
spaceery 27 U" 0 oad 
f > (F)G,t,) de 
= Tr) Wikre t E> 
space “on to) ° 
> ‘lk s iw(t-t) 
where 7 (fee y &r) =—— | 6 (rome dw . This expansion can 

0” 0 von 0) 7-0 


> 
be inverted as b(r,to) = ff 


a 
f> (ee y 2) T(t', rp) dtjun ae ether emi sea 
0 0 


unique one to one relation between and T as long as wW remains in 


the time domain D. Observe that 
Cee | :> @ a [ s (r)T(t',r,) at! 
ie = | Zz r r = ae 3 PEG § t 
C space Xo tt to) —00 is to) a 
= [ de" T(t" e) | Gy (r) dr 
00 ee space Tg (tty) T(t! -ty) 
= f dt! T(t!',z))8(t'-t) 
phis means o(t'=t) = f £> (eat yo) > (ict )@) de - This can be 
space 0 0 0 0 


verified by substitution if one remembers that only w > 0O is allowed in 
the distribution of energies. This in turm restricts the admissible 

eal 
functions T(t,r5) as functions of t , so the transformation above has 


the same properties as 6(t'-t) on the admissible time amplitude functions. 


Let us examine this transformation which we called "d(t'-t)" 


in some more detail to see exactly how it acts on time amplitude functions. 


yy 
pl a <i an 3a 


aah, fe ; hal 


“ as = 


_ OP 5. rae 7 Oe a | 


5 


-—— oa 
ey, 


nes notaeeqes ahiT ~« Wb 


—.% wie 

7 “an TIS SS 

se st etedT « “ab (gts F Gh a ‘)Q0 a} * Gay mee 

4 

. <7 a5 . : 

al sebasst v ee grot ex T bab Y asewied wanaaten ane ame supinw 
jactd svaeed0 


* - -_ < 2 
"46 gianna at | ” Ca ¢ amadQa 


| (oxo tap \~ : 
ihe 


“a Gh an 
lg? Sq?" 79" sana! 


wb Pe aes 


far'apd( a '3)T ‘ab _\ = 


} = (3-"3)? 
a7aqa a 


: ty 
od #ho stat 3b Ge Pe =" tg 4 OI a3) : 
at Bewolis sit 0 < w ¢iao seu oxaduenb: ono Tt votsuansedae wd a 

éidbeateba ad3 staivdest oxi? a} eff? .s0lgisee Qo meee iit: ie 

i 

eat avads qutisizotauerd 849 06 « J Lo senha? #8 (guat obs 
.enokisnm) sbuttiges outs ofdiantebs od] 00 (9-73) as - 
¥ . see 

: ost atu schmnxe- 


we =e 7 
Ss ! 


- _ 
o 


ta" aye" bollae sw. 
sin sm 
Pern wey vee TL 


53. 


Observe that 
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The Time Operator 
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We wish to have f-> (xr) as eigenfunctions of the time 
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Operator tT? with eigenvalue t . T> operates on position 
eat Lee 
Osa) 0° 0 
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Now let us look at the opposite order of operators, namely HtT> ier « (r) = 
1a yl bap OI 
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this last term separately. 
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The Zero Rest Mass Particle in One Dimension with Positive 


Momentum Distribution 


The Schrodinger equation corresponding to | E = cp, is 
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; 3 
if 5p (St) = -ifec SS (x,t) 2 i ry ees Sean The solution is obviously 


W = £(x-ct) . This can be expressed as a Fourier transform, 


hee Wee 
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The momentum wavefunction is $(k,t) = g(k) e for k>O. We have 
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g(k) = 0 for k<0O. 9 satisfies ih se chkd , k > O . The momentum 
distribution lock, t) |? = eee is non-zero only for k > O and remains 


constant in time. Also the position distribution Wate i = |€(x-ct) [7 ; 
moves with the speed c of light in the positive x direction and maintains 
its form. Since f is the Fourier transform of g , the uncertainty rela- 
tion holds for the spread Ax and Ak namely AxAk 2s, however both 


Ax and Ak are constant in time, there is no "time of minimum spread". 
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We assume f (and hence g) is normalized. Here Xy = <X>Q + This shows 
the centre of expectation in position space moves with the speed of light 


increasing with time. 
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Let us consider the normalized energy eigenfunctions os ¥ (x) 
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be obtained by Fourier transforming A(w, x,) as 


Less) = uf A(W, x, ) Seow dw . 
0 


1 ik(x_-ct) 
Prom this we have I(t,x,) = —— ve g(k) e dk = Vo W(x, st) : 
at Join 


7 ee) pet aon: 
va bas caged pe ieee 
estigat (unpayt = wb. =) Pi ' fo 


- td gt ee ee fe ay et 


Pe 
ap been ent G24, 
yt 


» 


“o (x) & oy = (2.20 


s] 


394 of xual! 10¥ Le aad _ am =o Jedd Tes 
( %e be « : 
— =~" G* Sas 
jag 02 WY @ot? 4 @ agiiszgeint ody ot 


Pe 
th BPO Fat xsl) By) | "3 (an 4 


w 
i | - 
: _ a «4 ya . 7 Cent of, peatm, ( eyladA oN = (a mey 


mao (,%43)T nataonvi stotiique sate saY . 02 uw i 0 — oe ine 
an 4% = guteiolens24 retro qd be ‘ 


ah “saad Fe ts 


ee ee ae Comte ‘ r 


gi 


T(t,x,) = yo £(x,~ct) - The relation between the time amplitude function 
and position wavefunction is T(t,x) = Vo W(x,t) . Let us evaluate the 
expected time of arrival at xX: This is given by 
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CHAPTER VI 


SCATTERING THEORY IN ONE DIMENSION 


We have considered the wave function of a free particle in one 
dimension. Now, let us consider a potential V(x) which is localized about 


x = 0. We assume V(x) = 0 for |x| Paes 


V(x) 


iar tha eh RAE 


Local Potential 


We then consider the Schrodinger equation 


2 2 
h 
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In scattering theory, a basic assumption is made about the 
position wavefunction wW(x,t) . We define the interaction probability 


of the particle with the potential at time t to be 


a 
p(t) = f lwx,t)|% ax . 
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We assume that (x,t) is such that p(t) tends to zero (or is essen- 
tially zero) for t outside of an interval which is the time interval of 


interaction with the potential. 
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Fieure 12 


Time Interaction Interval 


In otherwords, the particle does not react with the potential 
before and after scattering, and the solution for the wavefunction is the 


free particle solution. The momentum wavefunction before scattering is 


ink? (t-t ) 
= 0 
@ (kt) = $,(k,t)) ¢ 2 
and after scattering it is 
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b,(k,t) = oy (ksty) e ot 


oT is the wavefunction of the particle for t < Ty and b5 is the momen- 
tum wavefunction for t > To + Each has the free particle time dependence. 
7) is determined from >) by the momentum transformation kernel 
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$5 Ck, to) = ‘ oe a oe (1) 


for all t. . One can see that 6 


0 depends linearly on 6 


» since %, 
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and >, are respectively the asymptotic solutions for »? at t = —© and 


+00 respectively where $ satisfies the linear Schrodinger equation 
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which is linear, although the operator V(i ~ makes it difficult to 


solve. Alternatively, 6 is the Fourier transform of wW satisfying 
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Specifying initial conditions 6 for $ determines $ for all time, 
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and hence determines 7) - Moreover this dependence on > is linear, and 
so > depends on %) . This justifies the form (1) of the transformation 
shown above. Moreover we can see that the kernel Be Ccsk*) 3 like the 
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Fourier transform kernel -——e preserves norm on the transformed wave- 
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functions, since if we assume o, is normalized, then $5 must be normal- 


ized also (for >) a wave packet in momentum space). 
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and the position transformation kernel A, (x,x") by 
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In operator notation, equation (8) reads as 


The Delta Fiumction Potential 


Let us consider the simple potential V(x) = Vo O(x) , and 


determine the position and momentum transformation kernels. 
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Let us take the initial momentum wavefunction for large negative 
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Hence the final momentum wavefunction for large positive time t is 
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Potentials of Finite Width 


Let V(x) be a potential (usually positive) that is zero out- 


side the interval -a< x <a. The strength of the potential is 
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assume the potential is positive and does not change sign. 


The delta function potential gives a good approximation to the 
Momentum transformation if the deBroglie wavelength of the incoming parti- 
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We can do the same type of calculation in this case, assuming an initial 
and final momentum wavefunction of exactly the same form as in the case of 


the delta function potential. We can solve the time independent Schrédinger 


equation in this case, because we know that the time dependence is Bie 
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These are the explicit solutions for the coefficients of the 
position wavefunction. Because of their awkwardness, we shall simply use the 
symbols B,C,OD,E and A rather than the explicit expressions given in 


(10) above which are listed for reference. 
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By Fourier transforming we obtain the intial and final momentum wavefimction 
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It the. potentially EVKx)  is-localz~so- that Vix)= 0 for 


x S-a and x >a , the free particle equation is satisfied, and so the 


positive momentum solution takes the form 
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WwW Nas a behaviour dependins on V(x) for -a <x <a, but B and C 
can be determined for given V(x) by matching values and derivatives of 
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where C and B_ depend on Ko and B, (ks ko) is independent of t. 


Probabilities of Reflection and Transmission 


For incident particles (or wave packets) of momentum k, , the 
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In the case where A is real and positive we have 
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The transmission coefficient T can be measured as a function of ko 

In the case of the square potential, it provides a means for obtaining the 
parameters Eo and a. For the general potential, it gives us the absol- 
ute value of the functions C and B of ko - It does not give us any 
information about the phase, and consequently the momentum transformation 


kernel is not completely determined. We might ask if time delay measure- 


ments can provide this phase information. 
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Time Delay: hese Relation. tO. the Momentum Transformation Kernel. 


Suppose that the initial wave packet before interaction is 
described by the momentum wavefunction , (kt) with the free particle 


~ihk’ (t~t)) /2M 
time dependence >, (kt) = o, (ks ty) e » Let us Suppose that 


>) has a localized positive momentum distribution about k, > 0 , and o 
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wave packet interacts with the potential, is scattered, and the resulting 


free particle wave function is >, (k,t) . For a momentum transformation 


kernel given by 
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Figure 15 


Momentum Distribution 


After Scattering 


Determining the Potential V(x) From the Momentum Transformation Kernel. 


The WKB Approximation. Suppose that the initial energy level 


=e is large compared to the maximum height of the potential over the 


2 0) 
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range of the narrow momentum wavepacket centred at ky gar SG Suppose 


also that the potential is smooth and slowly varying over the range of one 


wavelength A = eh. 
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0 
V(x) 
— > 
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Figure 16 


Slowly Varying Potential 
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We use the WKB approximation to represent the solution (x) to the Schré- 
dinger equation in terms of the potential V{x) . We match this solution 
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From this we can express B and C in terms of the potential V(x) 
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Continuity of w and its derivatives at +a give us the equations 
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approximation, there is no reflection, and only transmission. In this 
case it is not the reflection, and transmission coefficients which provide 


us wich a direct tink to the potential V{x) , but rather the phase of the 
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transmission coefficient C that satisfies lc = oq ° 
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5 has the same momentum distribution as >) - Consequently measurements 
which are designed to determine y(k) as a function of k must involve 
position or time delay. Observe that (13) is a good approximation for high 
energy, and yet from it we can determine the potential V(x) which can be 
used in lower energy exact calculations to explicitly obtain the momentum 


transformation kernel. 


We note that both $5 and >, have positive momentum distribu- 
tions. Consequently, we can apply the free particle time operator as we 


have derived it in order to make time delay measurements. 
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so that $)(k,t) = $,(k,t) et ACk) 


Let us assume that ) has a very narrow momentum distribution 


about ko 2.0", and Jer ws -caleulate “<t> with respect to >) and then 
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Case 2: The Momentum Wavefunction 5 
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Now taking the infinite series for A(k) , differentiating and 


substituting we get 
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If we substitute Vf ae, into (14) we get the classical time 
delay result 
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The above calculations show that the classical time delay corresponds to 
the expected value of quantum mechanical time delay for high energy parti- 


cles and slowly varying potentials. 
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CHAPTER VII 


QUANTUM MECHANICAL TIME MEASUREMENT IN ONE DIMENSION 


The measurement of time of arrival in quantum mechanics is 
limited by fundamental difficulties inherent in time measurements. Unlike 
position, momentum, energy or angular momentum, time is an improper and 
not not a proper quantum mechanical observable, This means that an opera- 
tor associated with time measurements is not self-adjoint but only Hermitian. 
This means that time measurements cannot be made on every possible state of 
the system defined by a given Hamiltonian, but only on a restricted class of 
States defined by a subspace of the Hilbert space in question. The following 
free particle example illustrates this. 

Consider the free particle in one dimension. Classically, the 
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Hamiltonian is H = ae independent of q. Thus 7;—= =q and 7—=0= 
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In momentum space the time operator for a free particle in one dimension 


is 
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The unknown function (not determined here) turns out to be 
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f£(hk) = = where Xo is the position at which time measurement is taken. 


Thus the time operator is 


This can be used to make time measurements. The expected time of arrival 
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For narrow positive momentum distributions about Ko > 0 we can estimate 
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relation. 


The problem with T is the appearance of k and e in the 
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denominator in the terms - me ye and ea ae This imposes a restriction 


on the class of admissible momentum wavefunctions on which Ty cam operate. 
(e.g. it excludes momentum distributions concentrated about k= 0.) This 


means that To is not self-adjoint, but only Hermitian. Its domain of 


operation is not the complete Hilbert space. 


We can illustrate the statement that an operator conjugate to 
the Hamiltonian will not in general be self-adjoint but only Hermitian in 
the case of a Hamiltonian with discrete energy levels like that of the one 
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we know that ae ) is a Hermitian operator satisfying this relation. 


The Time Delay Operator: Its Domain 


By considering the momentum wavefunction at any specified time 


to we can apply the free particle time operator to =e eS AK Okt care 


to make time measurements at the position Xy ° 


As we have seen, the time 
Operator is Hermitian but not self-adjoint. If we construct the time 
operator for a potential and subtract from it the free particle time oper- 
ator, we obtain what is known as a time delay operator. If the potential 

is zero, for example, the time delay operator is the self-adjoint zero 
operator. It thus may happen that for some potentials, the time delay oper- 
ator is self-adjoint, or at least Hermitian, with a larger domain than the 
free particle time operator, for instance. On the other hand, the inter- 


section of time domains for the free particle and the given Hamiltonian 


may result in a smaller domain for time delay. 


Commutator of the Time Delay Operator with the Hamiltonian 
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ih = [H,t] = [Ho tH, »T,+T 4] = [Hy >T + Lez! + [H,T 4] 


ih + [H, tT, ] + (H,T4] 


This shows that [H,T,] = -[H 


Since Hy = V(x) isa function of x it would be useful to 
express the free particle time operator as an operator on the position wave- 


function. 


Operator Correspondence and Integration 


Our Hilbert space of functions of one real variable correspond- 
ing to either one dimensional position or momentum space, consists of 
functions (complex valued) defined on the range -~ < x < +0 with the 
property that when multiplied by arbitrary polynomials in x or differen- 
tiated any number of times the result is square integrable. This insures 
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element in the domain of -i De be a derivative of an element already in 


the Hilbert space. 


With this definition, we have the operator equivalences 
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The domain of tT consists of those functions for which the above expres- 


sion exists and defines a function in the Hilbert space. Observe that 
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This operates on the position wavefunction w(x,t) and corresponds to time 


measurements at the position Xo: lt.ds Hermitian but not self-adjoint. 
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Time Delay for Interaction with a Potential: 
Example: (The Delta Function Potential in One Dimension.) 


The potential V(x) = Vo 6(x) in one dimension has no classical 
analog. Classically this potential cannot be penetrated, and so measure- 
ment of time deiay on transmission is meaningless. Let us ask the question 
what is the time delay for penetration of this potential barrier in quantum 
mechanics? Let us choose our time domain in such a way that only positive 
momentum contributions appear for x > 0. Observe that approximation 


methods (like W.K.B. approximation) do not apply to this potential. 


The first thing to do is to obtain the energy eigenfunctions 
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By assumption there are only positive momentum contributions for x >? 0 
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This gives us the relation es Vo C(wW) = 21 [== Bcw) and so B(W) and 
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C(W) can be expressed in terms of A(W) as 
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It only remains to determine A(W) in magnitude and phase. The 
absolute value of A(W) is determined from the relation 
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zation is general for a local potential in one dimension. The phase of 


this result we have | Acw) |? = 


A(W) is determined in order to make the eigenfunctions positive and real 
at the point of arrival. We denote the complete phased, normalized eigen- 


function by WU (x) 2  dtvis eiven,by 
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¥ (x5) is a real function selected to make the phasing correct. We require 


w = (%) > Q . This can be divided into two cases, Xo > 0 sand Xo <a) 


The case Xo > 0 is the simplest one, and is of greatest interest, since 


we are measuring time delay on transmission through the delta function bar- 


Pret. LOY Xo > QO we have 


goietags «xb ‘abbey W GM | +a 


= Ve, i 
-biseron eiut ~ geisinonat 2928 . oe. - | enabine agit a : 
ae 


io Sealy sat AaM i pamepinte ps TOL3& 
inex bum ovisizoy naokyomvTanges: eda sole od! zabz0 mt wtae ae 
Ye eee ie oh inven ie ; 


~fspts bosilewrom -b9 beendg axbhiand ddd stansb oi sieviean q ait : 
re? : ® Ge ; et: 7 
of pavig as 3T ++ oat S pokao 
(evr 4 o 
5a, co ma 
; ——— # as Low, 
=e | - 1 - 
bas , O> x 
( x) ¥2 a 


A ; _ 7 : 7 7 ; 

ptlupes 2W =. 2597700 galasiq 3 elem 01 botoeloe notaomi Iset enh: i er 

eens ov? oont beblvib ed amo atdT - o<t — 
a ¥ aan 


an 


. 0? o* bas 0 < o* 
sonke ,32e93sIAl se93e0%— Yo Bt bnk ,900 jaslgaie ads ut o< 
a3gd cokiom? szish ofs dguetds nobsettenas? no yeiab oni a0 


pda 


2m mh 
iY (x5) = oo 0 (1 + 2 0 ? 
e = —— : 
nv. 
AL Sie 3 
20h 
and from this we can directly get ie = (x) for the case Xo > 0.8 01t. is 
b] 
0 
given by 
iV 
0 mh 
CL fa) 2) 
1 4 . n? 3u) a, Tae (xX) 
ve xX (x) = me Tw es : 
270 Jom 2 
mV) 
at 3 
2Wh 
4 iV 
Ae nm QO /mh — 
cpomsatd | samen GN a rie ae me 2m 
7 Jon 2 tw a WW ; dary A (xtx,) 
2 
ine [20 ca 
CEA) ee 
0 2wh 
£0f x00 | and by 
Ey ne 5} 
a 2hw tL (=x) 
y (x)= el i : tor Kye eu 
mVy) 
ik ae 3 
2wh 


In this time domain we are considering wave packets with a positive momen- 
tum distribution that are being scattered from a delta function potential. 
We can sketch the (absolute value squared of the) position and momentum 


wavefunctions before and after scattering. 
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[vox,t) | locke) |? 
delta function 
Incident potential 


Before scattering - position and momentum wavefunctions 


|wex,e) |? |oce,t) |? 
Reflected Transmitted 


After scattering - position and momentum wavefunctions 


Figure 317 


Delta Function Scattering 


Let us take a point Xo > 0 and look at the expected time of 
arrival of the transmitted wave packet for the delta function potential, and 
the expected time of arrival for the same incident wave packet at Xo ibis 
the potential V(x) had not existed. The difference between this two 
values for narrow momentum distribution about ko > 0 gives us the time 
delay in terms of ko . Let the wave packet be completely determined from 


the general position wavefunction 
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2 d hore: 
We assume |A@,x,) | (which is independent of Xo) is a narrow distribu- 
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tion about Wo So as The expected time of arrival at Xo ea0n 7s 


* 3 
<t> = [ A (w,x,) [-i a! A(W,x,)dW . Here we assume that (x,t) is 
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Recall that we wrote the energy eigenfunctions in the form 
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When we take superposition, the A (W) terms form the incident wave 
packet, B(W) the reflected wave packet and C(W) the transmitted wave 
packet. In the equation 
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for sufficiently large past time (negative t) only the incident wave packet 


of the eigenfunctions ve be (x) contributes. Hence if we let p, Gxt) 
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denote the free particle wave function asymptotic to (x,t) for t7*>-@, 
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Now let us ask what is the expected time of arrival for this wave packet 


for a free particle at the point X + Let ve ae denote the energy 
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eigenfunctions for a free particle in one dimension with positive momentum 
distribution phased at Xo ° Let AY W,x)) denote the energy distribution 


for vx, t) ;» then 
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The expected arrival time for the free particle described by p, (x,t) at 


the point x is given by 
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The time delay At for the interaction with the delta function potential 


is 


* ad rae fe) f 
= A (W,x5) [-4 3p! A(W,X_) dw - A (W,x5) [-i ae, A (Wx) dw 
0 0 
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for a narrow energy distribution |A@,x5) |? » about Wy: Expressing At 
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as the time delay for the interaction of a very narrow wave 


packet in momentum space distributed about k, > 0 with the delta function 
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potential V(x) = Vo 6(x) . If the range of the potential V(x) is much 
smaller than the wavelength a » this expression will be a good approxi- 
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mation in general. 
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General Relation Between Time Delay and the Phase of the 


Transmission Coefficient for the Momentum Transformation Kemel 


Consider a local potential V(x) in one dimension defined on 


the range -a <x <a. We have solutions w(x) to the Schrodinger equa- 
1 ll uote ~ik)x 
tion with the behavior (x) = ——e + B(kg) e for x < -a 
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for each ko >Q. For x > a we have only positive momentum contribu- 
tions indicating that the initial wave packet before scattering had only 


positive momenta in its makeup. 
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to note is that this normalization is independent of the local details of 
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[Acw) |? = an - Consequently, the normalized energy eigenfunction in 


our time domain is 
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The behavior of pO) on the interval -a < x < a is of no concern to 


us here. The only thing left to do is to phase the energy eigenfunctions 
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Now consider a general wavefunction wW(x,t) expanded in energy 
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Scattering interaction. Let ip, (x,t) be the free particle wave function 


whose wave packet coincides with that of (x,t) for sufficiently large 
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negative t . Only the incident waves in the energy eigenfunctions contri- 


bute to the actual wave function (x,t) before scattering, and so 


As in the case of the delta function potential, Gest) can be expanded 


in free particle energy eigenfunctions 
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From this we can derive a relation between the phase of the function C(k) 


and the time delay as follows. 
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At as a function of initial momentum, we can determine the phase of the 
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transmission coefficient C from A(ko) = 2 Ky At (ko) dk : 


Transmission Only 


The time delay At can be determined by an alternate method in 
the special case of transmission only (no reflection) in one dimensional 
scattering. In this case the momentum transformation kernel is 
BL (kyky) = ¥2m C(kp) 6(k-kp) and so $,(k,t) = Y2m C(k) $,(k,t) where 
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areival operator, namely tt -—"=——"+ et + a we can determine the 
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CHAPTER VIII 


NORMALIZATION IN THREE DIMENSIONS 


Let us consider the eigenfunction normalization equation in 


: * =e ie > 
three dimensions, namely | PCr v yr @) dr = d(w-w') . 
space 


As an example, consider the free particle in 3 dimensions. The 


>> 
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momentum eigenfunctions w>(r) ee are normalized in the 
k 3/2 
(2T) 
f Eg tk > > >> 
sense that | We) vey (x) dr = 6(k-k') . Now suppose we examine a 


space 
time domain and impose W normalization. If we take the simplest case 
imaginable (and one that is useful in scattering theory) namely the time 
AG ° . ae Bs 
domain of all negative momenta k directed along the z axis, k = -kz , 


k > O , we must normalize 


with an integration over all space. This is clearly impossible since 
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Wy 6t) is independent of x and y and these integrations will diverge. 
It is therefore necessary to consider solid angle distributions in momentum 


space. We write 


for the free particle in three dimensions. k is a unit vector specifying 
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the direction of momentum. We now substitute this expression for yy) 
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into the normalization integral | VOY @ dr = 6(w-w') to get 
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This is the normalization condition which must be imposed on the solid angle 
A > 
expansion coefficients A Ox) in the energy eigenfunction Wy) for the 


free particle in three dimensions 
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hold. The condition is necessary and sufficient. 
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Specifying A Ck) Satisfying the normalization coeffient for 
all 7 0 determines a time domain for the free particle in three dimen- 


sions. 


The direction of k can be specified by two spherical angles 
7 k) = 
kc and ae so that A (k) A(8, 29) Poatd 
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range of values of o. near o. = . For such narrow distributions, 
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where 6 satisfies 1+ cos 8, = Chen and [4,017 = Q for 


Now let us consider the actual value A ,(9,) . When we take the 
iy(w 3 01.) 
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Square root above there is a phase uncertainty of the form 
The W dependence is unimportant, since this directly relates to the 
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phasing of vee) which is to be determined later by phasing at point of 


arrival. The o. dependence is insignificant, since the distribution is 
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and it gives the amplitude of these waves in that region. We see for 

larger regions of plane wave behavior we have smaller ¢ and consequently 


less amplitude. 


mu) F - ° 
Bory ke= live and for momenta directed along the z axis in 


a positive direction we have 


sean aidan - BAS bl enn ot ontinaie b 


ae i &\ 4 es - - 

deg ai4* bel oss! #\4> lx) . xl x08 pied ah a 
s ‘ 

my 


(Syd Fo-abusi tga sz 9% 2 én Sud (bestepes Aw Ps 


bersi moos sw seafderqg odi enkolyes eidgT .c1ss 03 ylarokiow asa8 


~99nt Istunqa gateu « oidsiiav sm at * aot a: | 


jeavew sanity eft eaves Oy spins seftw 2ovo aet2iooge re 


140. 


ee meen) 2! lies 
Uae Sonn . 4 


for x , y , 2 satisfying the restrictions above. 


— 


ya 


met; ee 


i a 


; o ea 


CHAPTER IX 


EIGENFUNCTION PHASING IN TIME OF ARRIVAL 
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Let W(r,t) be expanded in evergy eigenfunctions as 
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and A(W, X59) on the time domain, we require that the phasing of the eigen- 
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function p (2) depends only on the eigenfunctions and not on the explicit 


wavefunction or energy amplitude function. 
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vector of maximum time variation. 


Suppose now we are considering scattering from a local potential 
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near = 0 in three dimensions, and the detector is a great distance 


from the scattering centre, and is off at an angle and out of the influence 
of incident waves that miss or do not interact with the potential. The 
detector region only contains spherical outgoing waves. Suppose we are 
measuring time of arrival at the detector. Let the detector be located at 
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CHAPTER X 


NOTES ON ORTHOGONALITY OF TIME EIGENFUNCTIONS 


Allcock! discusses the inherent problems in measuring time 


arrival in quantum mechanics. One of his observations is that it is 
sible to find a set of functions b, @) which are eigenfunctions of 


4 ve) 
time operacror §—i Ay and which are orthogonal in the sense that 
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[ PCO) PCW) dw = 6(t-t') . This is true in general, however if we 
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fact is, this integral does behave like a delta function (identity) trans- 


formation when it operates on the domain of time amplitude functions. 


T(t', x9) is a time amplitude function, then 
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Note the explicit use of the fact that A(W, Xp) = 0 “fom “o < 0.. =ihus 


1 Allcock, G.R. 1969: ANN. Phys. (N.Y.), 53,253. 
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if we restrict attention to the domain of time amplitude functions (this 
always happens when we work in a time domain) the required orthogonality 


holds. 
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The Space of Time Amplitude Functions and its Conjugate 
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for A(W) zero at W=0 and ©. Also yee) = [ (-iwA(W) ) 
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In order for this integral over W to behave like 6(t-t') we require 
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CHAPTER XI 


THE TIME DELAY OPERATOR 


In three dimensional space we consider the Hamiltonian for a 
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central potential H= - oM Ye + Vite Where tex |Z We note that 
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Classically, the time delay for a potential of range A is given 
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Gee esahe ally iat ap ade (tae aah aeand el ceeeMel at Gee edius 
2M Z 2 0 0 
2Mr 
we can write 
A 
dt = V2M 
. in 
m ae 
V 
The lower limits r_, and b of integration cause the denominators in 
m 


the respective integrands to be zero. 
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To make an operator out of dt we replace H and ‘ie by the 
respective operators. Since H and ie have common eigenfunctions 


we can see that ae are eigenfunctions of dt with eigenvalues 


A a A 
Nag = VOM [ Ee mH | 
(w, 2) 2 b(w, £) D 
hw - ACen - V(s) hw - Aten 
2Ms 2Ms 


The lower limits of integration in each case are selected to make the denom- 
inator zero. The eigenfunction equation dt = x comple 
oe q Yog m WR Yen pee) 


determines the operator dt since Wen are a complete set. 


Observe that H and ne are classical constants of the motion 
and quantum mechanically compatible observables. Suppose we make measure- 
ments of a state and find that the energy is always fw and the square 
of angular momentum is always R(R41)h- . Classically, the parameters b 
and Vo are specified and so dt is uniquely determined. For such a 


state any measurement of dt always yields the fixed value doe . «snus 


this state is also an eigenstate of the time delay operator and has eigen- 


value dg : 


Integrating the final term we can simplify the expression for 
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We observe that the time delay is simply a function of the two 


Z : ; 
constants H and L of the motion in the classical case, i.e. 


jeFSteG, Ee) 


A 
vm {= #8 _ ig 
Tin 12 
H - rue V(s) 
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b 12 
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Since Ho is not a constant, but a function of s we cannot make an 
operator substitution for H, and [i » even though there is a set of 


0 


eigenfunctions of Ho and Ee . Observe that 


Bt Wig = ECL = EChw, LLIB), o 


The time delay operator dt commute with H and - for a central 
potential V(r) . This follows from the fact that H , Le aid ote a1 1 


share the same eigenfunctions Yate . Moreover all the eigenvalues of 


these operators are degenerate in m. 
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Time Delay in Scattering 


We must enlarge on the concept of time delay if it is to be 
properly applied to the theory of scattering. In addition to our time 
delay operator dt we introduce the scattering angle operator O . 

s 


Recall that classically 


i} 
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i 2 
i 2 
wey Fs H - S - V(r) 
2Mr 
Tak: is selected so as to make the square root zero. Observe that Ee 


min 


is uniquely determined as a function of the parameters H and ie which 


are classical constants of the motion and quantum mechanically compatible 
observables. Thus the operator 0. satisfies the eigenvalue equation 


where 


p 


eaten = Tah 


2 dr 
~~ FE eceraTR | ee 


Oat 
- (52) | 2 
min a ne CAT h Shea 
2Mr 


Observe that the eigenvalues O52 are scattering angle values and are not 
directly related to the observed angular position of detection, since the 
direction of incidence in quantum mechanics may not be precisely known. 


Thus 0 is not an operator whose eigenvalues measure the angular compo- 
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nent of the point of detection. 


The operator dt can be used to measure time delay. If we 
measure <dt> on a particular wave packet for instance, we get a single 
average time delay value, but we do not obtain time delay expectation at 


a given scattering angle. 


The Time Delay Operator in One Dimension 


Let uS consider the time delay operator in one dimension. We 
Will find it is Hermitian but not always self-adjoint. In fact; it dis 
possible to determine for specific types of potential exactly what the 


domain of the time delay operator is. 
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Figure 18 
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Casey L: Nix): eerie Goes. 0, For |x| 2 a. The Hamiltonian 
ee 
ibs oe Soot V(x) has every eigenfunctions e.9) with H Ve:9) = 


dx —— 
. /2mw 
2 mw 4 m 7 h a 
hu Gx) ’ YO) er Cares Ph for x > a® , where c(k) is 
the transmission coefficient for the momentum transformation kernel. In 
(oe) 
x 
this way, YO) is normalized to satisfy | p 0) pd dx = 6(w-w') 
=e) 
FOr nw.) 0 eand pox) respresent a moneonergetic "wave packet" with 


energy hW which moves in the positive direction only for x>a. 


The p 6) » W > O form a complete normalized set of functions 
for expanding other functions in terms of in this time domain. A general 
wave packet which initially approaches the potential from negative x can 


be expanded as 


CO 
-iwt 
(x,t) = | Aw) ¥ (x) e 
0 
We can define a time delay operator on the Hilbert space spanned 
by Re:9) by specifying precisely how it operates on each b Ox) . We 
insist that By) be an eigenfunction of the operator ot with eigen- 


value Ay eythate 45 tp Cx)) = dey » where se is the classical 


time delay for a particle with intial energy fw to pass over the potential. 


Thus 
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Ma on ee ae 
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hep e.v. >0O be the maximum value achieved by the potential V(x) . Then 
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delay operator is that subspace of the Hilbert space spanned by Ww (x) 
W 


max 


oat It is a Hermitian operator, since its eigenvalues Ay 


rnoye (0) 


are; real but it is not self adjoint, since not all w(x) are in the 
u) 


domain of ot. 


Coseuc.: ViG@o) =O) -always, V(x)°="0 for |x| aaa. 


V(x) 


Figures 12 


Negative Potential 


Hyd is defined as in case 1. Instead of a time delay we leave a time 
advance in this case, since dw) defined above is negative. In this case 
at Cy (x)) = yy? defines dt on Y Ox) ope yeWuil | (hea Oh Vo) span 
a complete Hilbert space as long as the finite range attractive potential 
is not strong enough to allow the existence of one or more bound states 
(always finite in number) with discrete negative energy levels. If bound 
states exist, dt is only defined on that subspace of the Hilbert space 
generated by the positive energy eigenfunctions pox) for, W105.  CAny 


function in the domain of ot is orthogonal to every bound state. No 


time delay concept is associated with bound states. 
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If we take V(x) = Vo O(x) where Vo rat) ad can easily see 
that no time delay operator ot can be defined in this way since the 
delta function barrier cannot be penetrated classically and the integral 
for ro does not make sense. Thus the domain of dt is zero. If we 
take Vo <0, ot has a complete domain of functions b Cx) » (excluding 
the bound states since one exists), but the value of do is zero for each 
W, so dt is the zero operator. It becomes clear that our concept of 
what time delay actually is turns out to be less than ideal. The operators 
we have defined for 9t are known as classical time delay operators. In 
constructing them we have taken classical formulas for time delay in terms 


of energy in one dimension, and time delay in terms of energy and angular 


momentum in three dimensions, namely 


i dx 2a 


and assumed that these were definitions for the actual concept of what time 
delay means. In fact time delay in the sense we would like to think of it 
as, is not a simple function of energy or energy and angular momentum. What 
we have constructed is an operator for the quantity defined above in the 
integrals which is classical time delay. Therefore this operator dt 

should rightfully be called the classical time delay operator. We saw that 


in one dimension for a slowly varying local potential, with particle energy 
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much greater than the maximum potential height, the W.K.B. approximation 
told us that the expected quantum time delay was the classical time delay 
defined as in the above integral. First of all, we recall a calculation 
done earlier to determine the commutator of a time delay operator with 


the Hamiltonian H including kinetic energy and potential. 


The result was [HT] = am a - In our calculations with 
Ty = dt , we observe that and H have the smme eigenfunctions, and 
so [H,T,] = 0. In general, there is no reason to suspect why the poten- 
tial V(x) should commute with Te a free particle time of arrival 
operator. This indicates that the classical time delay operator ot is 
not the same as Ty the time delay concept, obtained by considering the 


difference of two time of arrival operators, one with the Hamiltonian 


including the potential and one for the free particle. 


A General Discussion of the Time Delay Operator 


Suppose we consider a typical scattering problem in which inter- 
action of the particle with the potential is limited to a finite time inter- 
val. Let the wavefunction be described by 

ng = > =iwt 
Crt) F switp We en (ee dw. 
0 a 
> 
W(r,t) remains in the full Hamiltonian time domain D _ spanned by the 


=< tJ ° 
functions ve = (r) as functions of Y for all time t . Initially 
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> > 
v(r,t) corresponded with the free wave packet pt) i p, (st) can be 


expanded in terms of the free energy eigenfunctions as 
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p, (r,t) remains in the free Hamiltonian time domain pf spanned by 


> 
ve is (r) for all time t . This time domain is a Hilbert subspace of 
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0 
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functions of the vector variable r. Since i ie (rt) is a multiple of 
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vy i (r) by a complex number of modulus 1 , the time domain pt is 
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c 
is and T> are operators on functions of position r. Observe 
PEAS es peyhs 
G0 0 
* aa 3> 
w (r,ty)T v(x, ty) r = the expected time of arrival 
space OmeG : : ae 
for a particle at point to in 
full Hamiltonian. This integral 
is independent of ty . 
* i 
w (%,t )T> ) (r,t Vd = the expected time of arrival 
1 On xt. 0 
space 0° -0 : ; ae : 
for a particle at point Ty in 
the free Hamiltonian described by 
(r,t) . This integral is inde- 
pendent of to - 
The time delay operator is T-> = cf for measuring time delay at 
ie SE EA sit 
0570 0-30 
=> 
To ° It operates on a wavefunction at time to which must be a time 
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before the scattering interaction. At that time v(r,t) = p, (r,t) and 
; : P P ; : f£ ; : 

this wave packet is contained in the intersection Dn D of time domains 


Since WeD and Y, € pt for all times. 
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at *z for full Hamiltonian 


paris time of arrival ad bee time of arrival 
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> 
at to for free Hamiltonian 


Time delay due to potential observed at ce : 


This time delay concept is different from the classical "impact parameter 
time delay" we discussed in the three dimensional case, since there we 
measured the time of arrival with and without potential at different points 
equidistant from the origin and separated by the scattering angle, while in 
this new time delay concept, the same a is used for both free particle 


and potential time of arrival. 


It is a good exercise in manipulating integrals to start with 
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CHAPTER XII 
THE DEVELOPMENT OF TIME MEASUREMENT IN 


3-DIMENSIONAL SCATTERING 


We already know the result for the time dependence of a free 
particle in Cartesian coordinates for two and three dimensions. These 


results are: 


os © 0 im(x-x)* 
¥Ot) = SaRCeEDy J, to) ex Gntee,) 03 


(for two dimensions x = (x) oX5))5 and 
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(for three dimensions x = (x) »X5 2X3). 


We shall consider the time dependence of free particles in other 
coordinate systems, namely 2-dimensional polar coordinates (which can be 
extended to three-dimensional cylindrical coordinates) and 3-dimensional 
spherical coordinates. We shall be interested in both position and 


momentum wavefunctions. 


Two Dimensional Polar Coordinates 
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Figure 20 Polar Coordinates 


In order to solve this equation for y we separate variables 
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Next of all, we use the Hankel transform representation of the 
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coordinates. 


Alternatively, let us consider the momentum wavefunction in polar 
coordinates. We denote it, by. 46(k,6.,t),..where; k = | & | the magnitude of 


the momentum vector, and £8 is an angle specifying its direction, so that 
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In the case of a free particle in two-dimensional polar coordinates 
we see the great advantage of working with the momentum wavefunction instead 


of the position wavefunction. 
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The simplest expression for the time dependence of v(p,¢,t) 
is found by transforming the two dimensional cartesian integral into 


polar coordinates. Thus we have 
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Separation of Variables in Spherical Coordinates 


For a spherical coordinate system, let us look at the free particle 
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Since spherical polar coordinates lend themselves well to a central potential 
V(r) which depends only on the distance from the origin, namely the 


positive scalar r, let us solve the energy eigenfunction equation 
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This last equation is known as the associated Legendre equation. It 

turns out that the only physically acceptable solutions occur when 
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integer for single valuedness. 
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at r= 0 is unimportant. The relation between the two solution sets is 
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P, (s) are the Legendre Polynomials, a form of the associated Legendre 
functions for m= 0. 

% is called the quantum number for total angular momentum, an 
m is the quantum number for z component of angular momentum. Ya (es 4) 


2 
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Recalling that f (8) = ) d, C)P  (cos6) we see that the second 
=o 


equation is a direct connection between the scattering amplitude and the 
function a» (w) which characterizes the information about the potential 


which can be determined from scattering. As usual, k = [2s , 


Although the asymptotic expressions are valid for large r 


we must keep in mind the restrictions Ix|,]y| << s/ Te and 
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lz | en | al . The dimensions of the incident wave packet are large 


compared to the range a of the scattering potential, but small 

P , , me Ae 
compared to the radius To associated with Xo the position of the 
director. Thus only the scattered outgoing component is measured by the 


detector. 


The absolute value of the scattering amplitude can be measured 
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experimentally by scattering probabilities, say aes | £, (8) | . We 
ask the question: Can time delay or time of arrival measurements 


determine the phase of £, (8) 2 


Scattering Phase Shifts 
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. B Cw) 

phase shift 6,(w) to satisfy the relation - tan 6&.(w) Eee oe es 
£ Q A, (w) 

(1) 


We can compare this with R yg = Cita, Cw) ho 


(ke)the”? (ker), and 
conclude that A, (w) = 2+a, (w) and By Cw) = ia) (w). Thus 


ia, (w) id, 
tan 5) =- ra, @) eo (that e = cos 6) + isin 5 = y lta, (w) 
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Thus we can express the scattering amplitude as 
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Consider a scattering interaction involving a particle described 


by a wavefunction 
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with a localized potential such that interaction with the potential is 
restricted to a finite interval of time, and only outgoing spherical 
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waves arrive at the detector located at ro: The incident waves are 
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not observed at the detector, and so the asymptotic form of Yo) at 


the detector is 
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Substituting for by @) and a, Cw) from their derived equivalences 
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define the phase mn, (8) of the scattering amplitude £09). Then 


i(kr.+%, (@)) 
= 

A(w,r,) = Ate © *k : 
Let us choose a particular form for A(w) which is narrow about W, and 


Square normalized. We take the dummy expression A(w) =/ 5(w-w) - 


YW, > 0. In practice, we have to take non-zero square normalized functions 


with a parameter that is non-zero, but selected to approach Y 6(u-w) 
as the parameter approaches zero. We then perform the required integrations 
for non-zero values of the parameter, and take the limit to zero after. 


In practice we can avoid these steps by the following formal operations. 
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We now use the relation A(w) = ¥5 (ww) to get 
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To obtain the time delay corresponding to 3t in the classical case we 
ask what is the expected time of arrival of a free particle wave packet 
with the same form as the given wave packet before scattering at a point 
of distance To from the origin along the z axis after passing through 
the origin. Before scattering, only the plane wave component contributes 
and so 
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x and y = To will not be in general). If we define ae for this 
case by z= to and x = y = 0 as the components for is and 


measure expected time of arrival we get (for the free particle) 
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Thus we compute the time delay ot a function of Ko and @ as 
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Observe that dt is independent of tg: 


In essence, we can determine the phase 2, (8) of the scattering 


amplitude from the time delay 3t(k,6). We have 
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m do h 
ot(k,6) = nie ak 1 68) and so n, (8) = a gt Ck, 6) kdk 
It appears that nh, 8) is not completely determined by this integral, 


but rather we can add an arbitrary function of 6, independent of k 

to hy, 8) and still have the same time delay dt. This arbitrary function 
can be determined by examining the behavior of the phase N,, £8) of the 
scattering amplitude £, (6) pi ee Op or k++. Under these 


d 2 , 
circumstances we can obtain m, €8) and a = | £,,(0) | from experimental 


measurements, and consequently we know £,.(@). It remains to relate 
£, (8) to the potential V(r) and consider the phase of £, (8) Tor 


limiting values of k. 
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Green's Function Solution to the Radial Schrodinger Equation 


Consider the equation 
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= C(r)sj, (ks) + D(r)sn, (ks) bOG,  <hjo eS > Or. 
G(r,0) = 0 => D(r) = 0. We choose a Green's function for outgoing 


spherical waves as s+ +e. Hence B(r) = 0. Thus 
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To determine A(r) and C(r) we have two more equations, namely 


3 af 


ht ae 


ahs orn go * a 
‘wv Pt y, 


as a932i2¥9" 


ve 
(a) that = coureae +2 


a 
MN te Shite: i 


t28)2 - Soy a paieny line! 
vnbsnups zononegomeeled oft conan Sa o+ oo" 
th . wind 
noliulor akan’ a'noesd sd7 epi! Ca)t = uti =a 
a) x 


1)3 = Ss eatisiadse (5, 17)9 ovsdlw woh [RN o1 m | 
+e | 
0s ~ Fy J 

_ "s 


apd? .anedsitewm yeubaaed suly ia-1)h = (Ce ’ 
' 


+ (a. 


7 - , 
(ot) aatare + cat Pita cada = (9309 al a 
a 


oe 163 

~O<8<?T WS Can) gatara i Cex) tets)9 of 
. ce 1 

i 
= ad 
gniegive yo! colyandl aieatS &.de0eis oW 0 = (200 =O 02 sinh 
‘polye 
wait .0 = (x) som et +e ae Seve aaah 


% = wt : 
33 cna" 40nd ~ (390 


> 
=i 


| 
v5 


| 
~s>4 
—_ - 


meet) eee» Ged) te) Oe 


a -_ 
G(r ue) = Ger) and Ag G(r,s) =]. 


From these equations we get 
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The factor in the denominator, r(j, (kr) kho 4 
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The Green's function solution is 
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This is the integral equation for the radial solution U (r). Furthermore 
W 


ae 
U (vr) eo 2r48 (kr) +o Cyan) Cee) for 2 > 4 and so fo ro 
w,2 Q Q Q , i‘ as 
a 
it 2 
ol diipan Sess er eee icone er rons 
Q 2 Q wr 
h oO 
Geos sib this integral, and&so G) (rs) = ~irskhe (ker) j, (ks) . Consequently 
we get 
a 
is 2mik : 
a, (w) = - 2 [. sj,(ks)V(s)U ,(s)ds 


The upper limit a can of course be extended to infinity. We assumed 
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The Born Approximation 


The approximate substitution U jg 6? = 2sj, (ks) into 
a, (w) = - 5 A sj, (ks) V(s)U , (s) ds 


is called the Born Approximation. It is valid in general for weak finite 


potentials and high initial energies. In this approximation, 
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Observe that a, (w) is imaginary negative for a positive potential 


2ikA, (k) 
V(s). Since at, (w) oye ee A, Ck) is negative real, and since 


foe) 


£, (8) = a A, (k)P, (cos6) F £8) is real in this approximation. For 

£,(6) non-zero we see that there is no time delay in the Born Approximation. 
The values of 6 which make £0) = 0 for fixed k are the angles 
identifying regions of destructive interference in the pattern of inter- 
ference fringes for <= | £0) |° - op any OG) vand) for which 

£,.(6) # 0 there is a non-zero probability of arrival, and time delay dt 
can be measured and will always turn out to be zero. If £6) = 0, 

the probability of arrival is zero and the time delay concept at this 


point for this energy is nonsense. This reflects our inability to 


determine the phase of £,.(8) At this point. 


Although ot = 0 in the Born approximation, this is not true 


of higher approximations or the exact solution in general. 


S-Wave (2 = 0) Scattering in the Born Approximation 
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The trick for inversion here is as follows. V(s) is defined only 
fon s+ Ua ex tendatherdetinition, of Vis) for = <0 by V(-s) = V(s) 


(even function). Then 
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in the Born approximation for the real scattering amplitude £,.(8). 
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Low Energy Scattering 


Suppose the potential V(r) has range a and ka << 1. Let 
solve the radial Schrodinger equation by Green's functions in this low 
energy limit. We have the equations 
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Since ka << 1, we can use the asymptotic forms 
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The W.K.B. Approximation for the Radial 


Schrodinger Equation 


Let us look for a semi-classical approximation relating the 


potential V(r) to the functions a) Cw) « The Schrodinger equation is 
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with a solution U g6®) satisfying U gS) = 2xj,(kr)+a,(u)rh 0” (kr) for 
r >a. We assume the potential V(r) “has the following very special 
properties: 
(Gee ee ees bounleds at )s<9V(0)| “forsall 7 = 05) -¥(0) 
finite. 
(2) V(r) defined on 0 <r <a is strictly monotone decreasing 
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This form is sufficiently general for present purposes. Under some 
conditions, one or more of the above restrictions may be relaxed. 
The approximation we will consider is one of large k (high 
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the solution is of the form (boundary condition at r= 0 imposed) 
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the variable of integration to k we have 
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We can compare the above equation with the classical time delay 
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the first direct connection between classical time delay and approximate 
scattering phase shifts. Note that a zero phase shift 5, (k) = 0 means 
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Vere ) = 0 andso V(r) % ad(r) which is the zero potential with 


an infinite peak at r= 0. 


The Second Born Approximation 
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U yg 6X) £Gre Pre ae namely U ig 6X) = 253 (er)+a,(w) ra? (kr). From the 


integral equation 
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wid errncis "| a 


~ Gal 2 
ab (a), UCerWKed) tz i =. c a ee 
+S 


css) ft ; : - 
-s > 3 went, tere) ew || pret aes _ 


we 
oy we -_ 
> & : 


alfa o¢mi (24), e% > (a)y,8 ainstjedve ae r¢ a _ 


ee a 


it) : sf. - : 
(aa) da tae . t 7 
ebtepy {isd ,f)°e eee a - 2 é ad a | 


etapa pegs) ptt aha 4. er. L 


ric booase pila mk (ayo dlebaneton oft 0 24 
osnk ovols muvig ea (a),.U arene 


CJ _ 
7 : 7 - nm 7 
a _ i 7 a iy Oo _ Oe 


a 


a, (wy = - 2mik [* |, (kr) V 
g ne soos (x) U(x) dr 
@) 
[e@) fexeie 
Gab Atle | 12 2 
a, Cw) laa, Te £ [j,(kr)] V(r) dr 
fh fe) 
DEEP STIS 1 
— | ae rj, (lx hy (ker) VC) | 8E5, (ks) ]?V(s) 8] 
h fe) fe) 
2a2ana a 
- = | dr MERC OMATON s*ny (ks) 5, (ks) V(s) ds] 
h fe) r 


If we write the latter integral in the plane and interchange the order of 
integration, we see it is the same as the preceding one and so ot, Cw) can 


be given (in the second Born approximation) by 
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We assume that the Born approximation holds, and so the first 
t erm above, although small, is much larger than the second term. 
Observe that ne (kr) = j, (kr) + in, (kr). We are interested in obtaining 
the approximate value of the small imaginary part of a, Cw) and the much 
smaller real part of a, (w) - Since the first Born approximation gives 
us a reasonable value for the imagingary part of a, (w) , we ignore the 
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imaginary contribution of in, (kr) and simply replace (kr) by j, (kr). 
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This means we can explicitly integrate the second order term and we have 


the approxiamtion, 
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The condition for validity of the Born approximation is 
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this first term is small compared to the first term but of the same order 
as the much smaller second term. Also observe that we can reformulate the 
condition of validity of the approximation for large k by using the 
asymptotic form (for large arguments) of j, (kr). The condition looks 
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The W.K.B. Approximation in General 
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to the argument of vi which has units of energy. We can change the 


variable of integration to Ww to get 


W J 
66 / a) = — v¥2mh i (eau) yee (hw')dw' - ka 


wW q 
= lim [- v2mh | Games (hw')dw' - kal 
a a v(a) 
h 
a d -1 ules : 
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This expression can be inverted as before, and the result is 
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W.K.B. Approximation for Larger 2% > 0 


For %&> 0, the Schrodinger equation for U Ae?) becomes 
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that R, (¥) = Vy (hu) where V is the modified potential 
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Bacibiatory wand for rs < Ry (w) U(r) degenerates to zero as r>QO. 
Let us assume we have a real function U(r) satisfying the differential 
equation above and the boundary condition at 0. We will ignore the 
normalization factor on U(r). The standard W.K.B. approximation across 


Classical turning points (with the connection formula) gives us the form 
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U(r) = A, (k) rj, (kr) a By (k) rq, (kr) for r > O we observe that 
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- tan 6, = —~ defines the scattering phase shift §& .(k) and more- 
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Phase Error in Approximation 
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1 
kr - VECRt1) z We know 6, (k) = 0 for V(r) = 0 and consequently 
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This gives us the fundamental form for the scattering phase shift 


in the W.K.B. approximation 
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Notice that this expression also holds for 2 = 0 , and 


6, (k) as given by this expression coincides with the same value obtained 
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earlier. A crude estimate for 6, (k) can be made assuming -v," (hw") 


contributes the largest values for w' near zero. Doing this we get 
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specified, V(r) is determined, and so therefore is 5) (k) for each 3; 
If measured values of 5) Ck) are not self-consistent it means that the 
Scattering interaction involves more than just a simple potential. The 


differential scattering cross section « = | £,.¢0) |? measures the 


magnitude of £,.(8) and the time delay can be used to determine its phase. 
With £09) explicitly known, the scattering phase shifts can be determined 


from the formula 


1 id) (k) 
£,.(9) = es (2L+1)e sin 5) (k)P, (cos 6) 


or78 


from which we get 


16, (k) 
e sin 6, Ck) = 


Ny 


T 
[ £ (8) P (cos®) sin 0d. 


Such a simple method for obtaining the phase shifts 5, Ck) for 


scattering would not be possible if only | £,,(6) | was known and the phase 
in (8) — £,, (8) 
of the scattering amplitude was undetermined. For e = £6) we 


we have (taking real and imaginary parts above) 


T 
Bk : cata 
cos 6, (k) sin 6, (k) aie i. | £,. (8) | cos h, (O)P (cos 6)sin 6d@ 5 


TT 
sin’ 6, (k) -* [ |£,,(0)| sin ,(0)P,(cos®)sin ode . 


Since (8) >+>nt as k-2>o (in the Born approxiamtion) we have 2, (8) 


can be found explicitly from the time delay as 
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Optical Theorem and Total Scattering Cross Section 
l oo 18 Ck) 

Peom@ ce, (6% a 1 (2R+1)e sin §(k) P,(cos8) we 
can evaluate the total scattering cross section o= | $3 dg’ where 
(0 ae 2 : 
aaa |£ (9) | . the result is 

o = a ) (2241) sin’, (k) 
Q 
k L=0 
l co id) (k) 
Also f, (0) == ) (22+1)e sin 6 (k) since P.(1) = 1. Hence 
k copie g Q 
“ Im(£, (0)) = 9 . _thas result is known as the optical theorem, and it gives 


an expression for the total scattering cross section o as a function of 


k in terms of the scattering amplitude. 


The Separable Radial Potential 


In the radial Schrodinger equation ee + Ge - yu ) = 
Ww 9 
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—> V(r)U @) , and in the corresponding integral equation 


Ung (e) = 2x5, (ier) +5 | G, (r,s) V(s)U,,, (s)ds 
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radial potential. For such a potential, the three dimensional Schrodinger 


equation is 


2 co 
SP bat h 2 > A 
tn 2b Cre ea SPs wrt) + | = Vir,s)(sr,t)ds 
Oo 
where t= rr and [x | = Fre 


We are especially interested in those potentials V(r,s8) which 


can be factored as V(r,s) = V,@)v,(s). Such a potential is said to be 


separable. 
From the integral equation for U gs) we have (taking Va 
and @V(E3s) = VCS) Vi(s)), 
Ue) = 2r4 (kr) , 2m [" Gears) ViGe)ds ‘ V(q)U .(q)dq 
Wh ig zZ J R a WL 
fh fe) fo) 
2m [- . 
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h E fo) 
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The equation above is an explicit solution for U gh) if N is known. 
N can be obtained by multiplying the equation for U sg? by Ve); 
integrating over r from zero to infinity and solving form N. (Since 
V(s) is assumed to be known.) Let us consider the separable delta 
function potential” Virvs)= Vir)VCes), = ve 6(r-a)6(s-a) so 

we Vo 6(s-—a). In this case WN = Vo U ig @)> We can evaluate Ug 


iti two. cases Tf <a and tr > a to get 


:. : Zmirk., Ci) 

Ug = 2rj, (kr) a jy (kr) VU (aah, (ka)V, tome oan 
7 ; 2mirk (1) : 
= 2rj, (kr) - —~ ho (kr)V QU, (a) aj, (ka)V, for ESA. 


2 
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If we set r =a we can solve for U jg (® to get 
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For the more general separable potential we can still obtain 


a (w) = - Busy s j.(ks)V(s)ds where N= ViiqgvVuleataddg . This 
Q n- 15 g ‘ Wk 


transformation is of interest because it can be inverted. The Spherical 


Bessel transform or Hankel transform is defined by the equations 


g(r) = iB | j, (ker) £(k) kak and © £(k) = /= | j, (kr) g(e)r“dr : 
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Thus we can solve for V(s) in terms of a Cu) to obtain 


Aree a 2 
pens, : hk. okdk 
WiGa) oy Sa k eae Neco eng 
@) == | 10a 
O R 
Thus for the separable potential V(r,s) = V(r)V(s) we can obtain the 
potential factor V(r) form a(w) or from the phase shifts 6 (k) 
: se ue 
216 Ck) 
where e = et Cu) « As in the case of a local potential, the 
phase shifts for a separable potential satisfies a consistency relationship. 
TE 6 (x) is specified tor ali k, alt phase shifts 6 Ck) are known. 
In the case of the separable potential, unlike the local potential, the 
expression for V(r) is exact and not an approximation. In particular 
V(r) can be explicitly determined from the s-wave phase shift as 


2 co w 
_ ih sin (kr) hk’. 


The usefulness is limited because the dependence of N on k 
is not known in general unless V_ is known. 


If we consider V(r,s) = V, @)V, Gs) where vy and V5 are not 


oO 


necessarily the same we have N Ce) = | V,(q)U (adda, with 
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2mik - 
a» Cw) = - i. Ny (k) | Ss j, (ks)V, (s)ds , 
h ) 
2 0 2 
> ith ? hk, _kdk 
Wee es | jy Cerda, 5) N, Ck) 


Suppose that a, Cw) = 0 for all mw and some fixed 2. Then 


Ug = 2rj, (kr) since N(k) = 0 and this would imply 


co 


N, Ck) = I V,(q)U Ca) da eae | Vo @)rj, (kr )dr 


jo) 
i} 


This means V,(@) = 0 by the Hankel Transform. Hence for a 
non-zero separable potential for every & there exists a k for which 
5k) # 0. The only transparent separable potential is the zero 


potential. 


Discussion of Transparent Potentials 


Let us see if there exist any potentials V(r,s) so that if we 
solve the equation 
2 tee) co 
Ue) = 2rjeGkr) te 16 G (rea) VelqU.(q)daide 
Wh vs per Z Q " Wk 
we get the simple solution ee? = 2rj, (kr) so that a, (w) = 0, 


For this to happen we must have 


0 = | CG (@,s)I | V(s,q)qi, (ka)dalds 
Oo (6) : 


at least holding for large r (in the region of phase shift measurement). 
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Since such r>s,where s is in the range of the potential, we can 


replace G, (r,s) by G, (r,s) = irskh)"? (kr)4, (ks). Hence we require 


oO 


Oa= | 23 ,(ks)1| V(s,q) qj, (kq)dq]ds 
@) (@) 


Since this holds for all k by the Hankel transform inversion we have 
0 = | V(s,q) qj, (kq)dq. Again by the Hankel transform inversion 0 = V(s,q). 
fe) 


Thus there are no non-zero radially symmetric potentials with 5, Ck) = 0 


for all, ky and any given £. 
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in the local potential case, we simply can replace V(s)U 6S) by an 
Ww 


operator on Ue in the non-local case and so 


yy 25 y@y 2) 
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Q we Q oll le Bie 
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fe) fe) fe) re) 


But the integral involving can be inverted as usual and hence 
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for all w > 0. Completeness of the radial scattering states (provided 


(emer C2) 


there are no bound states) implies V is the zero operator. 


Insofar as the scattering states are concerned, yh) = v6), 


An important corollary of this result is that all information 
one needs to know about a radial potential, local or non-local is 
contained in the s wave scattering phase shift 6 (kK) - This corresponds 
to the classical case for a radial potential in 3 dimensions where we 
found that the potential was determined by knowing the time delay 
9t(0,v_) for impact parameter b = 0. The s wave phase shift is 
given by 
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and this is the only quantity one actually needs to measure experimentally 
‘ do 2 : 
Knowing he [e203)) only is not good enough} the phase of £, (6) is 


essential. 


The Forward Scattering Problem 


Let us consider the case where £ (8) = 0 except in a narrow 
neighborhood of 6= 0. We ask the question: Is it possible to construct 
a potential V(r,s) such! that £, (9) Satisfies this condition for all 


k > 0? The total scattering cross section is o(k) = | le ,.¢9) |° ag, 


and the solution for £, (8) satisfying this condition is 
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£, (8) = / 2. a0 neces) e K By the optical theoren, 


o(k) 0 Im(£, (0)) => Te = Ove sin fact ru, is proportional to the 


narrow width of the "delta function" in some sense. Hence Im £, (6) = 0. 
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for all 06 meaning 6, (k) = 0. Hence £, (8) = () «and se.) aUk)) = 0. 


This means the potential V(r,s) = O by our theorems. 
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Three Dimensional Non-Radial and Non-Local Potential 


In this case the integral equation for the energy eigenfunctions 


+ 
y 62) in normalized form is 
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From the expression for £.(2) above, it is evident that if 
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for alle ck oand o This follows from the fact that the three dimensional 
Fourier transform operator is unitary with the full Hilbert space as its 
domain. The Fourier transform does not send non-zero functions to zero. 
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orthogonal complement of D. In this case EG) = 0 for all. kc and all 
rt. This potential operator V which is non-local and not radially symmetric, 


is said to be transparent with respect to the time domain D. It is clear 


that asymptotically at large r (where all experimental measurements are 


made) UG is the same as the free particle wavefunction 


£ meet PCy ened kz 
vj =e A k e 


It is also clear thet no scattering cross section will be found 
<o - 0, and time delay measurement cannot be made since ott) = 0 and 
cannot be phased. Since all measurements in any experiment are made 
outside the potential region, and outside the non-asymptotic region of 
the wavefunction, this potential cannot be detected at all. However, all 
is not lost. This non-local potential cannot be spherically symmetric 
as we saw, so we simply have to rotate the apparatus, and send incident 
particles in at a different angle. This amounts to changing the time 
domain from D to D'. If D' is in the non-zero portion of the domain 
of the operator V, scattering will occur, and ii) will not be 
identically zero. To get full information about V we may have to send 
incident particles in from many different angles. If V maps every 


time domain D to zero, V is identically zero on all scattering states. 


Any potential can be detected by appropriate scattering. 
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which is the free particle energy eigenfunction. Thus for a transparent 
potential with respect to a time domain D, the wavefunction is completely 
uninfluenced by the potential even in the region close to the potential. 


The general wavefunction with time dependence is expressed as 
eo 
> > -j 
wrt) = A(w)p (x) en dw 
fo) 


in terms of its energy distribution. For any (non-transparent) potential 


the initial free wavefunction corresponding to v@,t) before scattering 
> - —j 
p, @,t) = | ACw) vi eye’ dw, has the same energy amplitude A(w) if 
Oo 


phased at a point ahead of the scattering region. The magnitudes are 


always equal regardless of phasing point. 


Suppose the operator V describes a potential that cannot be 
detected by scattering, that is, V is non-zero only when operating on 


its own bound states. We can prove that V = 0. 


Since vy @) ="0O™ for all. do andvalt ime domainss. DD tor 
the incident wave packet, the scattering amplitude is identically zero 
> f > E \ : f> 
and Y = Y 2) in each time domain D. Hence Vy = 0 for 
all w>O in each time domain D. But the free energy eigenfunctions 
see Se 
£ > ; . : ik+r 
yW (rt) for all time domains, include the exponentials e for 
Ww 


all k and span the complete Hilbert space. Thus V=0. This justifies 
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the statement that every non-zero potential can be detected by scattering. 


Recall that the integral equation for the scattering states is 


> ess m 
v2) = Vy = eeaoy 


and the scattering amplitude is CO 
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CHAPTER XIIT 


THE NARROW WAVE PACKET APPROXIMATION IN ONE DIMENSION 


Consider a wave packet in one dimension approaching a potential. 
Suppose that the potential is slowly varying and the wave packet has nar- 
row width which remains narrow for all times at which the (position) wave 


packet interacts with the potential. 
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[WGx, ty) |, ty>t, 
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Motion of 
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bigure 22 


Narrow Wave Packet 


We assume that the potential V(x) is essentially constant over the non- 

zero spread Ax of the wavefunction (x,t) for any fixed t such that 

<x> is in the range of the potential. Im particular this means we have 
12 


transmission only, and no reflection. The Schro dinger equation 
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3 hao 


ih WS 05) Sealer (x,t) + V(x) W(x,t) at any fixed time t can be 
ie 2m xe 


approximated by replacing V(x) by the constant Vics) = p(t) 2 The 
constant potential (in x) is accurate over the narrow region of non-zero 


~ in position space. Thus we get the new equation 


2S 
Perro lee 


This time separation of variables gives equations which can be explicitly 


solved in general, unlike the Schrodinger equation. Take wW(x,t) = X(x)T(t) 


pete b): e Mees Cc) eit 
and get ih T(t) p(t) = hw and oS este oe = fw . Thus 
F = fae 
ikx mh : : 
X(x) =e and T(t) = exp[- 7 p(s)ds - iwt] . The general approxi- 


mate solution for (x,t) is a superposition of separable solutions, so 


that 


; ; F te 
W(x,t) = i‘ NOs en aa Ne =| (sy ds Wide 


a 


In terms of the momentum wavefunction 


Ke 


iki hen (Et. ey eet: 
(x,t) = ae $(k,ty) 2m exp [- =| p(s)ds] dk 
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Recall that p(t) depends on the particular wave packet chosen, and that 
Ax the position spread must be small compared to the range of x over 
which significant variations in V(x) occur for thosetimes t at which 


Sx? is in the range of the potential. Also note that p(t) = V(sx>) 5 
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: 1 ikx 
since. Ux.) = ib Ost) eo . dk wetcan compare and get 
V2T7 2 
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-ih Loe (t-t 


2m 0? 


A ec 
o(k,t) = (k, to) S exp[- i | p(s) ds] 
iE 


0 


Since lock, t) |? is independent of t , we expect this to be a 
high energy approximation. From the time dependence of the momentum wave- 
function we can get the momentum transformation kernel. We assume p(s) 
is non zero only for s in some finite interval, and so for to Be ad OS 
t*> +° #4=we have 

2 


beige 
-ih ae (t-t)) i 
o(k,t) = o(k,tp) e exp [- Lf p(s) ds} 


ane 
Sor ee ra 
Poreany t. 4 to + -o we have o, (k",t") = o(k' to) e 
2 
-ih _ (tt) 
Similanly £0r any t& , © +7 we have >, (k, t") = o(k,t) e 
We observe that $, (k,t") = [" BL 1(k,k") o, Ck", t") dk' and so 
2 “02 
-if A (ct *-t) -ih A (t'-t5) 
o(k,t) e = = f Bi 1(k,k") o(k',to) e dk' . Sub- 
Z 
-ih £ (t-ty) ; 
stitute $o(k,t) = o(k, to) e exp[- if” p(s) ds] to get 
2 
gua = oy i ere 
(k, tp) e exp [- a p(s)ds] = [ BL (kk )>(k sto) 
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Ber (k,k") = v2 C(k) 6(k-k') . Substitute above and simplify to get 


C(k) = exp i= a P,.(s) ds) =. 


This expression gives the function C(k) which determines the 
momentum transformation kernel in the narrow wave packet approximation in 
one dimension. The approximation has lead to an odd result. We expect 
C(k) to be a function of k but independent of the form of the wave 
packet. Actually we find C(k) does not depend on k but depends on the 
function p(s) which is p(s) = ASE which depends on the average 
momentum of the wave packet that determines the rate that <x> moves 
through the potential. The subscript k on P,. (Ss) refers to this 


average momentum. 


Let us now decide what type of function to choose for the initial 
momentum wave function o, (kt) which has free particle time dependence. 


We choose the minimum uncertainty product wave packet we derived earlier, 


namely 
21h (eecks)? ike 
RG 2 iin ©1150-5231 Gt ao 2 
1 4 (Ak) 

where 

ty is the time of minimum spread of position wavefunction, 

<k> is the mean (average) or expected momentum (constant in time), 

<x? is the expected position at time t) 


Ak is the momentum spread (a constant in time). 
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The mean (average or expected) energy 


z 2 2 
1S* Sheltie (<k> (hk) | The 
2m 
position*spread “Ax “at time’ t*"is 
Terenas: \s 
2s 2 1 
(Ox) ~ as 2 
4 (Ak) 2m 
Set t= to in the expression above for % (ty =) and wse 
2 
-ih *- (t-ty) , ft 
o(k,t) = >, (ky t ee exp = al p(s) ds] to get 
—0O 
ar aes 2 rae rt 
2 CCL Ne MR Oe See eee a oe er OL 
uf 4 (Ak) co 
4 fc 
= (kt) expl- ¢ [. p(s) ds] 
Thus we have the momentum wavefunction $(k,t) determined. 
Including the dependence of p(s) on average momentum we have 
ihe 
o(k,t) = $, (k,t) exp [- a | P,.¢s) ds ] Now let us evaluate en 
<1. cs on the function 
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‘3 h t ay ge 2 UngeP 
-<of +2 do(t-t,) + f —— exp esis) (| — = (s) ds] dk 
1 -© VT nk 2 (Ak) ao 


We can integrate by parts in this last expression to get 


2D, 
ES) exp tS) | 


fo ooh 
<x>) = <x> + m <k>(t-t,) ee i { i P,, (s) ds ]dk 2 


1 co /2T t(Ak)> 2 (Ak) 
Thus 
q h (k-<k>) Gicccl>)s 
ae <X> mee icone peamepe aece ilbs eae rte P(t) dk 
-0©/2nm h(Ak) 2 (Ak) 
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28 eee 
= <k> + [ |o, (kt) | fs ( yk oaks 


This gives the deviation of <x>) from the free particle case 
due to the influence of the potential. Recall the expression we derived 


for expected time delay for a local potential in one dimension, namely 
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Let us examine the conditions needed for the narrow wave approx- 


imation to be good. Let us suppose the potential V(x) has range 2a 


Cive. “@Vtx)-=20 for | x| Say Let sb = the maximum value of | V(x) | 
|v" (x) | 
and suppose that ae. is = ko for some ko ail: Ko finite (i.e. take 
max 
aye ere he 9) | a 
ko iia max | = | ). Then ko» eae. ae are the three limiting prop-— 


erties of the potential V(x) we wish to consider. 


Now let us look at the parameters specifying >, (kt) and see 


how these relate to the conditions of the approximation. 


(1) Momentum spread Ak >O. 


(2) Average momentum (<k>) > 0 . We require (<k>) >> (Ak) since 
the incoming wave packet is assumed to have only positive momenta 


contributing in the superposition. 


(3) Time of minimum spread ty - We set ty = Q for convenience. 


i 
oO 


(4) Expected position at t, . We take <x> This means the 


i 
free particle achieves its minimum spread at what would normally be the 
centre of the potential for the potential absent. This is reasonable. We 


want Ax to be as small as possible in the interaction region. 


In estimating the condition on the potential V(x) for the Narrow 
wave packet approximation to be good, we ignore the effects of the potential 
itself on the free particle wavefunction p, (kt) . We assume that the free 
particle wavefunction does very little spreading, and Ax remains at the 
minimum value (approximately) permitted by uncertainty throughout the passage 


of the region -a < x <a. We also assume that Ax << a and 
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; . i 
Ky Ax <<l1 , i.e. Ax << min {a , —-} and remains so. Recall 


2 2 2 
1B CAk)“(t-t,) 


us 
(Ax) | = ay oe ae The spread is a minimum at t = ty 0%. 
4(Ak) 2m 
; = z _ ma 
The wave packet centre is at aa =a at time ty = ners The 
2 
h (Ak) ?e 1 9 <k> 

controlled spreading condition is Sa) he ee ee GANS GS | Goer oe 

2m 2 (Ak) * 

Ko 
The condition Ax << a means (Ak)a >> 1. Also mee 1 comes from 
ky Ax << 1 , and of course we had (<k>) >>JNle¥) Thus 
a 1 1 1 1 
See eet ae —— —_ ‘ 
Eas SNF min ta , kp and => << Ak ° It is easy to see that the 
: : ? 1 
parameters a, ko - ire of the potential itself must satisfy ko > = 
by applying |[V‘(x)| oy ee ky to a triangular potential with height 
V at x =0 and satisfying V(-a) = V(a) = 0. Thus min{a , SS Ae 
max ky Ko 
Thus oe Kk. << Akh< es <k> . Without Ak we have wee ce 
a 0 oa a a 0 

Se << (<k>) . Al] these inequalities hold if we merely assume ky<< “e - 


ik 
Since =< ko for any potential, we derive (<k>) >> = and hence 


JS <x (<k>) . Thus the only condition we need to impose on <k> in 
: ‘ : <k> 
order that the narrow wave packet approximation be good is ko 4 tay a 


In addition to this high energy condition, we must also have the momentum 
: : <k> Sa 2 
spread Ak satisfying Ko << Ak < aii The condition (<k>) >> ko a 


is both necessary and sufficient for the narrow wave packet approximation 


to. .be-gend. SLEY (<k>)Le> Ke a simply take Ak = Se and construct 
the free packet p, (kt) with t, = O and Sa ees 0 according to the 
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formula for the minimum uncertainty momentum wavefunction. In order to 


insure that the variationsin potential over the range of Ax are small 


Cons f< sti) armom s >> xh notsthane | we 
eodt . ah <<. (<a?) nna 


asl esans If >> 


i a 
i*s” Ca.) 7 ape 


2« A yialsse taum ifseti Intstasjog onaz Yo aa . vf « 5 isaad 


silt 3904 soe Oo YRSo BP JT 


tfytort dtiv Istinsseq rsluganks3 ao e = 2 |@'¥| antulags ah 

° oe = be . Slain ait . O = (CoV ie~)7- apr vieizace bos Ox 35 a 

s a ae 
3 gis 4und ow JA aworleiW ais >> | oe 

3 Box seuees Ylotem ow 14 biad aslstisepant seeds (fA tees > S 

somen bie - <¢ (<a>) seylreb ow ,letansz0q Yas 702 al 7a & sont’ . 


al <i do seuget o8 bean 34 aokstbnos qino sa euil ) (<i>). —e 
; gi = beng od cotventxetq¢s tedeeq evaw worsen 943 sod 020 
euamomcs oft svad cely Jeum ov »wotsihnos Yas9aR dgid atd?.o3 cotanhinet, | 
Pa 3 << (<a>) cok sibuos od? . Bs > AA >> go? gatylatise 4A mie 
notdsalxoxqgs fsdoaq svav worsET od 207 Jostoliive bas seine dood & 
saus3e002: bas =\ = 4h adaa vigete ne Aci) Te -bacg od 0 0 


nee a ss 
r yr _ ; 
madcap ha ann ete, ee 2 6 | 


Kiem ore xh 40 span enn 5 al 2 
: . 9 ae 


‘ - wy cer 


sj 63 gatisosds 0 = 


bp) 


219. 


2 2 
compared to the incident energy we require beske 7aV » that is we do 
2m a (max 
Se ae 2 ee 
h <k> h <k> 
A ae Fe Case, ee ee ° 
not have mn Nees In fact we require om o we in order to 
construct the classical formula for P,, ft) namely | 
P(t) = Viy,(t)) 
where 
~] m dy 
YosC9 o> : 
cs a 0 2nV (y) 
pp oe 
hk 
dp, Ct) 
In this case x = 0 when t = 0. Let us evaluate annie for constant 
A i. . = ! == = i 
Upseinis is given by ig Vv" Cx) ak where x Y,.00) and t is 
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dp, (t) 
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3 op, (t) 
We will consider evaluating [s, P(t) J dt. Having evaluated ewotae 


: ; dt 
at constant t , we multiply it by as dx where ee is taken at constant 


op, (t) 
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Now integrate both sides. The two terms on the right hand side can each 
be integrated by parts. The limits on x may be taken as -a to a or 
—- to © , In either case V(x) vanishes at both limits and so there are 


no boundary terms in the partial integrations. Hence 
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Let us look again more closely at the approximation conditions. 
The spreading condition on Ax in the range of the potential was derived 
by replacing the potential itself by a free particle. We have no idea 


what effect this has on the approximating conditions. It turns out that 
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the condition cal Pa: will have to be replaced by placa V . 
2m my — max 2m max 


The first hint of this appears in the time delay expression 


te a V(x) dx 
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differs in the aa term. The best way to show the breakdown of the narrow 
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wave packet approximation for low energies on ~—oBeaT aes isto 


consider total energy. Recall that the momentum wavefunction in the narrow 


wave packet approximation is 
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i 
(k,t) = o, (kt) expi[— 2 P,.¢s) ds] . The expected value of total 
(oe) 
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This gives us <V(x)>, = 0 so clearly this is not much of a potential. 
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To give a consistent approximation, Cole << a For a narrow 
wave packet BAA San ~ Vox) = Pep (t) . The approximation is guaranteed 
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is the differential equation satisfied by $(k,t) in the narrow wave 
packet approximation. If we multiply by 6*(k,t) and integrate over k 


we get 
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E is constant, and f |o(k,t) |? P(t) qi Kehoe 
—0O 


Like the Born approximation, this is a high energy weak potential 
approximation and not semi-classical. For a more complicated semi-classi- 


cal treatment of narrow wave packets see Lebedeff. 
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APPENDIX I 


THE HERMITE EXPANSION FUNCTIONS 
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The Hermite expansion functions are po) = : H (x)e ~* We : 
n 
2 alt 
If F represents the Fourier transform operator defined by (Ff£)(x) = 
1 ikx = 
— f(k)é die hens Fye= dep < Bad dey = Cio WP with 
n n n n 
2 *=—3 
Yo 
vy 
respect to the basis » F is an infinite dimensional diagonal 
Vo 
matrix with i” in the nth diagonal position. (The initial position 


corresponds to n=0 and not n=1.) Clearly F is a unitary trans- 


formation. If we represent the function f(x) = ) f Ve:9) by the 
n=0 
fy fo 
fy if, 
column vector £ » the Fourier transform of f is oe . We repre- 
2 2 
e if, 


sent linear operators by infinite dimensional matrices and functions by 


column vectors in this system. 


Some of the properties of po are the following. 
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the formal expansion 
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Summation and Integration Conventions. 


There are two types of variable used in this analysis, discrete 
and continuous. A discrete variable takes on only the values 051 ,2,3,4,... 
and a continuous variable takes on any real value. If a symbol representing 
a discrete or continuous variable appears on both sides of an equation, it 
is understood that the expression of equality holds for all allowable values 


of the variable. If a discrete variable appears on one side of an equation 
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only it is understood that a sum from 0 to © is taken on the term 
containing this variable. If a continuous variable appears on one side 
of an equation only, it is understood that the term containing this 


variable is integrated from -© to ™ with respect to that symbol. 


As examples of the above conventions we have 
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The Hermite polynomials HO) which are used to define Vy OX) 


satisfy the following conditions. 
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He) - 2t H(t) + 2n H(t) =) 


on or > (n-1) 


Z = = 
wee y (<1) 9n ZY Ae 2r 
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The differential equation for p Gx) can be obtained by 
applying (x + Sx = =) and is given by 
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> + (2nt1 - x yw = 0 
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The delta function O(x-x9) is given by 6(x-x5) = do WU, 


From this we can show that 6 (x-X) is an eigenfunction of the operator x 


with eigenvalue Xo > that is x0 (X-Xp) = XS CeX,) . 


. Of eee 


ae As) 


g250%4 (USMerr tr id 
ey ** * or ae 


fe @.* . ar : PP 4) 
Bey, eb - 3 
, 1m 
. Wes) = Ye Jee wop * | 
< ‘= - @-s? ai eae 


qd Bentetdo ad ase, oe 707 nakinups tstaosast32b, sdf 
yi novig el bas &- gt” 


a= gd canat eg Aad 


on he L = igre ot aaah a peer 


Proof: 
x0 (x-x 9) = ) Wx) xh, GD 
n=0 
‘e a n n+l 
= 2 v Gt /® eg ORD ee ese Me (x) 


= ty ee Os age ae 1) <2 ob w | (XQ) ee (x) 


ime yil 


= 1 Ya FY + a v1) [2 9,00 


. 5 fort 2 
7 Z Yi+1 Xo) + [2 We BOS. 


(oe) 
=) [xGh (9) 1b (x) = x95 (xx) 
n=0 
Notice that the sums above are divergent as they stand and each term is 
understood to represent the kernel k(x, Xp) of some operator on a class 


of functions which allows the sums to converge. 


As a special case of the above result we get xd(x) = 0. This 
can be used to determine properties of the derivative of the delta function. 


Aside from the important result 
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where [n] denotes nth derivative, we can obtain by differentiating 
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result to obtain information about a class of 


functions we call delspan or general. First of all observe an important 


analogy about raising and 


lowering operators, orthogonal, analytic and 


general or delspan representations of a function (or generalized function). 


Three classes of function (operator) expansions are given below. 
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As an additional result we have the exponential expansions 
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The first two we know already, while the last one is proved later. 


Fourier Transform of a General Function 


Let 2(x)as ) b RTE AS be general with expansion coefficients 
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The function space S is the largest one which allows all sums and integrals 
to converge when a given function is explcitly operated on from S. The 
integrals and sums representing the operator themselves do not necessarily 


converge. 
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v2 
Thus 
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function is analytic, and the Fourier transform of an analytic function is 


general. The orthogonal functions 
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are both analytic and general. 


Expanding an Operator 
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Hence g = Ken te with summation convention of m. In fact using the 
Summation and integration conventions, we can write the above equations 
coe Ray) eee, hs ity) == 5 DC) aot) Cy) 

g@) = g UG). 6% 8) Ue), eg 5k £, and g(x) = k(x,y) f(y) - 
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The equation oi a fa is the matrix representation of g(x) = k(x,y) 
E(y) 4 The matrix elements are oe = k(x,y) wy Gs) py) . The general 


operator equation is ki = 2. 
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Suppose an operator A is expanded in the form 
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Can every operator be expressed in the form shown above for A ? 


The simplest separable operator has the kernel 
(n,m) x ‘ 
E (x; Vv) = ae) en) for some fixed n,. and m. The matrix elements 


(n,m) (as n)ex : ; 
of E (x,y) are Ey. = ae yoae - The operator k with matrix ee 


can be expanded in terms of the separable operators as k = ae ee . As 


a kernel, k(x,y) = k _ (Pm) 


sie (x,y) with summation of course. Hence if 


the separable operators can be expanded in the form shown for A , any 
operator can. The operator A can be written in the equivalent form with 


different coefficients b as follows. 
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Clearly, if A is expressible as in (*) it can be expanded also 


in the form 
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Also linear combinations of operators of the form (*) are also of the 
form (*). Similarly linear combinations of (Tt) type operators are (T) 


type operators. Moreover composition of (Tt) type operators are (tT) type 
operators. The operators x and < are, both of type (4) and (71). 1£ 
we can prove that compositions of (*) type operators are (*) type operators 
then we know all (+t) operators are (*) operators since every (1) operator 
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both operators on the right hand side are of type (*). Hence operators of 


type (*) and type (+t) are equivalent. 
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Any term in the sum for which n->m is zero, since in that case 
Cae =)" pbs) = 0. We can recognize this by the appearance of a nega- 
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Hence every operator is of sehe forme(*) or.(7) . 


Usefulness of the Operator Representation 
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entiable or not. 0 and T agree on the domain of VD” , but T is an 
extension of aS in a very natural sense. Since the domain of pk 
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includes all polynomials, the functions, e PX for all real DP. and 
We(x) sfor n= 0;1,2,3,.0-| wnich are examples, of ‘complete sets, we reter 
n 
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The Operator ee 5 


One type of linear operator is function multiplication where 


f(x) is transformed to y(x) f(x) . Another very special type of opera- 
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; ikx 
Expanding e as a General Function. 


Recall a(x) = 1 b, 6!"l(x) where pb = CY 
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a nl ee dx . Thus 
.n P eri 
Ee = and x) . This gives e “= 25 J a Seon i Meee 
; n=0 ~~" 


we can expand the kernel of the Fourier transform operator as 
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This is to be compared with 
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It is important to remember that (4) is an operator equation and not a 
converging sum of numbers. In order to verify completeness of the general 
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In turn, it is sufficient to verify this equation holds when 


the functions a(k) and b(x) are basis functions from a complete set, 


since the equation is linear. We choose a(k) = pads and b(x) = ig ike 


The left hand side is 
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This proves (4) and hence the general function expansion is a complete 


expansion. 
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But the coefficient in the expansion of a general function is always given 


in terms of the integral of the expanded function. Using this formula we 


get, 
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Fractional Differentiation 
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The Fundamental Coefficients 


Assume summation and integration conventions apply. Define 
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APPENDIX II 


THE FRACTIONAL FOURIER TRANSFORM 


Recall that from Appendix I some properties of the Hermite 


Expansion Functions are: 
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This series is divergent, but we are to understand that the 
functions f to be transformed by this kernel are such that the coeffi- 
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APPENDIX III 


PARTICLES UNDER THE INFLUENCE OF A POTENTIAL (HARMONIC OSCILLATOR) 


We shall consider particles in one dimension under the influence 
of a potential, and will solve the problem completely for. one special case 


(the harmonic oscillator). 
Consider the equation 
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where the potential V(x) is independent of t . This is the general 


Schrodinger equation which determines the time dependence of the position 
wavefunction W(x,t) for a particle in one dimension under the influence of 
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We attempt to obtain specific solutions by separation of variables, 
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Depending on the nature of the potential V_ solutions may be found for 

various values of E. The values of E for which solutions exist may 

be discrete or continuous. For the free particle we saw that a continuous 
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Harmonic Oscillator 
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This is the position wavefunction for the harmonic oscillator. The coeffi- 
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Expectation values to be considered are 


To do this:we must recall some.properties of the Hermite Expansion Func- 
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The Momentum Wavefunction for the Harmonic Oscillator 


The momentum wavefunction is determined as 
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The Time Shift Operator for the Position Wavefunction with 


the Harmonic Oscillator Potential 


We have seen that 
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Here, k is the fractional Fourier transform operator. We can substitute 
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This proves that the harmonic oscillator wavefunction is periodic with 


2 
period a Caste Ses cy ee Ca 44h - This can also be seen from 


the expansion of the time shift operator, namely 
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Time Shift Kernel for the Momentum Wavefunction 


Similarly we have a time shift kernel oe ) (kk) for the 
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momentum wavefunction @(k,t) . It is given by 
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